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Stability problems (buckling, critical speeds, etc.) arise in numerous elastic 
systems. The typical situation is as follows: Under sufficiently small loads, a 
configuration g of stable equilibrium exists in the immediate vicinity of the 
equilibrium configuration g) of the unloaded system. For certain (critical) 
loads, however, g either becomes unstable or removed so far from gy that the 
system is endangered. 

During the last decades numerous solutions have been worked out and 
compiled (see for example [1]*), [2], [3] or parts of [4]). Yet the different methods 
of solution, though their interrelation is but loose and by no means evident, 
never have been compared nor critically analyzed. Thus it could happen that 
lately a number of apparently simple problems showed up for which the differ- 
ent methods yielded different, in some cases even utterly improbable, results 
[5], [6], [7]. The reason for these discrepancies was found in the fact that the 
systems concerned were nonconservative [5], [8], and it was possible by im- 
proving the basic assumptions to obtain trustworthy results [9], [10], [11]. 

The incident showed the desirability of an extensive analysis of stability 
problems and of the methods for their solution. Such an analysis, based on a 
classification of mechanical systems, is the subject of this paper. In order to 
avoid unnecessary complications, it will be restricted on the whole to holo- 
nomic and scleronomic systems with differential equations that can be linear- 
ized?). 


I. Mechanical Systems 


In order to proceed to a systematic classification of mechanical systems, we 
start by circumscribing some simple concepts which usually are not defined 
with the rigour or the generality desirable for our purpose [12]. 


1. Work 
The work done by a force F is usually defined by the integral : 


W =| Far, (1.1) 
i 


taken over the path C of its point of application. This definition, based upon 
the concept of a stationary field force F(r) disengaged from the body on which 
it acts, is lacking in precision. 

In the first place, many forces are instationary, i. e. they are functions of 
the type F(r, ¢); and other forces may depend on the state of motion, especially 
on the velocities v of their points of application, and thus be of the form 


1) Numbers in square brackets refer to the Bibliography, page 49. ‘ 
*) The author is indebted to Prof. R. GRAMMEL, Prof. M. Rauscuer, Prof. E. P. Popov, Mr. 
M. Beck and Mr, Cu. Wenrti for valuable help and numerous excellent suggestions. 
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Fir, v). In both cases (1.1) is useless so long as the motion is unknown.’ When- 
ever r(¢) is prescribed, however, the quantities y, F and the power P= Fv 
become functions of the time ¢. Thus, W is better defined by the time integral 


w= [Pa=[Fod. (1.2) 


Example 7.7, If a particle m (Figure 1.1) is free to slide along a rod rotating 
with variable angular velocity w(t), the radial component of the centrifugal force, 
V,= my w(t), is instationary. 


ee I. A 
eT 
AE Ue 
w/t) r 
ay 4 
Figure 1.1 Figure 1.2 Figure 1.3 Figure 1.4 


Particle sliding on a Projectile. Rolling wheel. Band brake. 
rotating rod. : 


Example 1.2. The air resistance acting on a projectile (Figure 1.2) depends 
on its velocity; it may be described approximately by R= —dvv. 

Example 1.3. The friction force F acting on the particle of Figure 1.1 depends 
on its relative velocity; if N denotes the normal pressure, F = —p N v/v. 

In the second place, many forces (such as reactions) are not defined apart 
from the body on which they act. If the dimensions of this body are finite, the 
displacement (or velocity) of the material point on which Facts may be different 
from the displacement (or velocity) of the geometric point of its application. 
In any application of the concept of work (energy theorem, principle of virtual 
work, equations of LAGRANGE) dr actually stands for the displacement of the 
material point; thus, (1.2) is to be completed by the prescription that o relates 
to the material point of application of F. 

Example 1.4. In the case of a rolling wheel (Figure 1. 3) the material point B 
is instantaneously at rest (while the geometric’ point B moves). Thus, the reac- 


tions N, F do no work. 
Example 1.5. In the case of a band brake (Figure 1.4) the total amount of 


work done by the (apparently resting) forces 4, B is negative, since their mate- 
rial points of application move. 


According to (1.2) the work done by F over an infinitesimal displacement is 


dw=Pd=Fvd=Fdr. (1.3) 


2. Conservative Forces 


The concept of conservative forces is introduced in order to single out 
those systems which comply with the theorem of conservation of energy. 
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Usually, a force is called conservative if its work W is independent of the path 
connecting two generic positions of its point of application. This definition, 
too, arises from the idea of a stationary field force. 

According to 1, systems subjected exclusively to stationary field forces are 
extremely scarce. Thus, a more appropriate definition, valid for any (and even 
for reactive) forces, is based on the condition that the work done by F over 
any admissible displacement uniquely depends on the initial and the final con- 
figuration of the system. By an admissible displacement, in this connection, is 
meant a real (in contrast to a virtual) finite displacement of the system, com- 
patible with the constraints. 

Let 41, %, ... denote the generalized coordinates. Suppose the system to be 
moved from a generic configuration 4, g,,... to a standard configuration 
dio = 920 = --. = 0. Provided that F is conservative, its work is a single-valued 
function W,..9 = V (4, %, .--) of the coordinates, the so-called potential energy 
of the initial configuration. Hence, the work done over an infinitesimal displace- 
ment is 


dW =—dV = 3". Pidqey Pld ds a) ==> (ee 


On the other hand, the radius vector of the point of application of F is 
r(71 42, ---, t). In rheonomic systems (characterized by moving constraints) it 
depends explicitly on ¢, while in scleronomic systems (with steady constraints) 
it is a function of g,, J, ... only. For an infinitesimal displacement 


or Or | 
dr = ) 04, dq; a OL at. , (2.2) 
Hence, (1.3) yields 
0 
dW = SF 5 qe FS dt. (2.3) 


Comparing (2.3) with (2.1), we find that, provided F is conservative, 


Or 


i 


Or : 
Fa e='0, F i ae (¢=1,2,..,) (2.4) 


The first equation (2.4) implies that, apart from rare and insignificant ex- 
ceptions, conservative forces do not appear in rheonomic systems. 


Example 2.1. The particle of Figure 2.1, moving along a uniformly rotating 
rod, is rheonomic, since its position with respect to a fixed coordinate frame 
depends on y and ¢. The vector 0r/0¢ is normal to the rod and of magnitude wy. 
Hence, any force not directed along the axis of the rod is nonconservative, for 
instance the normal pressure N, the work of which increases the kinetic energy 
of m as it moves away from O. 
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We now are justified in restricting ourselves in the further discussion to 
scleronomic systems. Here, since P; and r depend on qj, q», ... alone, the remain- 
ing equations (2.4) imply that, in general, only forces depending uniquely on 
G1, Ja, -.. are conservative. 

In any system, two kinds of forces can be distinguished (see [13], p. 12). 
Active forces are given a priori as functions of ¢, 4), Ja ---) G1» Jo, ---» Whereas 
veactive forces are unknown beforehand, being obtained together with the 
motion only by integration of its differential equations. 


Or 
y CL Ne yy Yh C=2m oY 
Za in 
0 wt bas 
Figure 2.1 Figure 2.2 
Particle sliding on a rotating rod. Coriolis force. 


V, in Figure 1.1 and R in Figure 1.2 are active forces; Fin Figures 1.1 and 1.3, 
N in Figures 1.3 and 2.1 are reactive forces or reactions. 


Some reactive forces never do work in scleronomic systems; thus, they are 
conservative. The work done by any other reactive force in a real motion is 
negative. Forces of this type are called dissipative forces (see [14], p. 257) ; they 
are nonconservative. Thus, reactive forces are either workless and hence con- 
servative or dissipative and hence nonconservative. ; 

Normal pressures and static friction forces (see [13], p. 37) are conservative 
reactions; kinetic friction forces and rolling friction couples (see [15], vol. I, p. 133) 
are dissipative reactions. 

As for the active forces, there follows from (2.4) the rule that those depend- 
ing on the time or on the generalized velocities are nonconservative. 

If the particle m of Example 1.1 is displaced over dy at a time #,, and over 
—_dy at a later time f,, the total amount of work done by the centrifugal force V, 


generally does not vanish. 
The work done by R in Example 1.2 is negative for any actual motion; thus, 


air resistance is a dissipative, though active, force. 

There is an important exception from the last rule. Some active forces, 
called gyroscopic forces, are always perpendicular to the velocities of their points 
of application. Though they depend on the generalized velocities, they do no 
work in actual motions; in consequence, they are conservative. 


Example 2.2. In a rotating coordinate frame, the Coriolis force C=2mvu,x @ 


(Figure 2.2) is always normal to the relative velocity v,. 
Other examples are the gyroscopic moment of a symmetric gyroscope (see 
[16], p. 345) and part of the Lorentz force in an electromagnetic field. 
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On the other hand, active forces depending only on the configuration of the 
system are not necessarily conservative. Since the coordinates 4, Jz, -.. are not 
determined entirely by r, some of them are not even stationary field forces. 
In order to distinguish forces which, though functions of qj, 9, ... alone, are 
nonconservative, from those depending on qj, G, ... or ¢, they will be called 
circulatory forces, the notation being borrowed from a similar case in hydro- 
dynamics. 

Example 2.3. A constant force F, applied at the point A to a rigid boa 
(Figure 2.3), is a conservative stationary field force. 


oa ce 


¢ \ 
/ ie 
! \ 
A = 5 
AAy | 
tt Je 3 
\ 
\ 
\ 
\ 
Figure 2.3 Figure 2. 4 Figure 2.5 
Constant force acting Force rigidly connected with a body. Constant couple acting 
on a rigid body. on arigid body. 


Example 2.4. The force F of Figure 2.4, rigidly connected with a body and 
therefore constant for an observer taking part in its motion, is not entirely deter- 
mined by r; thus, F is not a field force. When the body is moved to a second 
position by a simple translation, F does a certain amount of work. If this trans-_ 
lation is preceded and followed by equal but opposite rotations, the work of F- 
is generally different; hence, F is circulatory. 

Example 2.5. Let M be the constant moment vector of a couple (Figure 2. 5) 
applied at O to a rigid body which is free to rotate about the axis z of M. Provided 
that the angle of rotation does not exceed 2 2, M is conservative. For angles in 
excess of 2 a, however, or if, instead of the axis z, the point O alone is fixed, M 
is circulatory. Two successive rotations through 2 about the axes x and y result 
in the same final position as a rotation of equal magnitude about the axis z. For 
the simple rotation the work done by M is different from zero, whereas, for the 
compound motion, it vanishes. 


Circulatory forces are more common than is usually supposed. Since they 
do work on bodies moving perpetually through the same positions, they play 


an important role in power-transmitting devices (such as shafts, se tae crank 
mechanisms and pulleys). 


3. Conservative Systems 


A mechanical system is called conservative whenever all (active and reac- 
tive) forces acting on it are conservative. This definition is practically useless 
without the precise definition of a conservative force given in 2. 
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From the analysis of 2 follows 

Theorem 3. A scleronomic system ts conservative if 

(a) tts reactive forces, in actual displacements, do no work and if its active 
forces are made up exclusively of 

(6) noncirculatory forces depending on the configuration of the system only, 

(c) gyroscopic forces. 

On the other hand, nonconservative systems include 

(a) rheonomic constraints, 

(0) dissipative forces (active or reactive), 

(c) circulatory forces, 

(d) instationary forces 

apart from more general (but practically less important) forces depending on 
the velocities. 

In a conservative system the work done by all forces while the system is 
moved from a generic to a standard configuration is a single-valued function 
V(qy) 4 ---) of the initial configuration. It is called the total potential energy 
and can be obtained by addition of the potential energies of every single force. 
For an infinitesimal displacement the total work, according to (2.1), is 


OV 


57: @=1,2,..) G2) 


aW = J) Pdq,, Pp = 


The principle of virtual work, applied to the actual motion, yields the energy 
theorem dT /dt = P, where T denotes the total kinetic energy. By integration 
and introduction of V the theorem of conservation of energy, 


T+V=E = const (322) 


is obtained, valid for any conservative system, whether gyroscopic forces be 
present or not. 


Example 3.1. A particle (Figure 3. 1) moving under the sole action of its weight 
is conservative. This proves right even if the motion is referred to a uniformly 


M 


Figure 3.1 Figure 3.2 Figure 3.3 


Relative motion. Double pendulum, loaded Spherical pendulum, loaded 
tangentially. by a constant couple. 
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rotating coordinate frame, since, beside the weight mg, the centrifugal force 
(m w? x, m wm? y), as well as the (gyroscopic) Coriolis force (2m w 7, —2m w #), 
is conservative. 

If w depended on the time, all forces except mg would be instationary. If 
air resistance were accounted for, the particle would be acted upon by a dissipa- 
tive force. In both cases (3.2) would no longer hold. 

Example 3.2. The double pendulum with elastic restoring moments (Figure 3. 2) 
may be regarded as a simplified model of an elastic rod built in at one end. 
Unloaded it represents a conservative system. However, the load F, directed 
along the axis of the lower rod, is circulatory (Example 2.4); thus, (3.2) is not 
valid for the loaded pendulum. 

Example 3.3. A spherical pendulum (Figure 3.3), subjected to a constant 
couple of moment M, is another nonconservative system, since M (Example 2.5) 
is circulatory. 


4. Classification of Mechanical Systems 


According to Theorem 3 the only conservative forces appearing in sclero- 
nomic systems are workless reactions, gyroscopic forces, and noncirculatory 
forces depending only on the configuration of the system. With regard to 
stability (see II), the presence or absence of gyroscopic forces is important. 
Therefore, a conservative system shall be called gyroscopic or nongyroscopic 
depending on the presence or absence of gyroscopic forces. 

On the other hand, the more important forces rendering a scleronomic 
system nonconservative are of the dissipative, the circulatory and the insta- 
tionary types. Transferring these notions, too, from the forces to the systems, 
we obtain the classification of Table 4. 


Table 4 
Classification of mechanical systems 
(with examples of stability problems) 


Mechanicel systems 


Saenens snwerenery 
Pp 


( 


Euler's Critical speeds Buckling under Buckling by Buckling by 
buckling cases | of unloaded influence of @ tangential 2 pulsating 
shafts damping load load 
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Obviously, Table 4, in its fifth category, includes even rheonomic systems. 
The only systems excluded are those containing forces which depend on the 
generalized velocities, belonging, however, neither to the gyroscopic nor to 
the dissipative types. With respect to elastic stability, they are of little impor- 
tance. For further investigations, systems with negative damping perhaps 
would be the most interesting ones. 

It is clear that, although a system is either conservative or not, it may 
belong to more than one of the last five categories at the same time. 


Example 3.1, forinstance, with variable w and air resistance taken into account, 
is gyroscopic, dissipative and instationary. 


In order to avoid unnecessary complications, we shall confine ourselves, as 
a rule, to pure categories, tolerating, however, in all of the five categories, the 
presence of nongyroscopic conservative forces, since, in elastic systems, they 
always appear (in the form of internal forces). 


Thus, by the term «purely gyroscopic», for instance, we exclude the presence 
of nonconservative forces. 


Table 4 contains five typical stability problems. That they all belong to 
different categories, 1s easily verified by comparison with Examples 2.3, 2.2, 
ieope 4, and Isl. 


5. The Equations of Lagrange 


A mechanical system is called holonomic if the differentials of its coordi- 
nates (and thus its generalized velocities) are independent. 

Nonholonomic systems are relatively scarce. A well known example is the 
wheel of figure 5.1, rolling with horizontal axis on a horizontal plane. Whereas 


its coordinates ¥, y, x, pare independent, their differen- 
tials are connected by the nonintegrable relations 


d%#=vcosadm, dy=rsinadp, 
representing two nonholonomic constraints. 

In a holonomic system the increment 6q; of a single 
coordinate g,; represents a displacement compatible with 
the constraints. The principle of virtual work, applied 
in turn to these displacements, yields the equations of 


Figure Dail LAGRANGE, 
Wheel rolling on a 
horizontal plane. d of OL _ Os en a aes (5.1) 
dt 04; 04g; 


the so-called generalized forces Q, being the coefficients of the increments 69; 


in the expression 


Ww = SQ; 69: (5.2) 


ZAMP IV/7 
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for the virtual work. The differential equations (5.1) are valid for any holo- 
nomic system, whether conservative or not. 

Any virtual work is evaluated under the assumption that, for an instant, 
the actual motion of the system is interrupted and replaced by the virtual 
displacement. The forces acting on the system, however, are supposed to re- 
tain the values corresponding to the actual motion. Hence, the coefficients Q; 
in (5.2) are not necessarily identical with the P;’s of (3.1). In (3.1) every term 
belongs to the same (actual) motion; in (5.2), however, the Q; belong to the 
actual motion, while all of the other terms are associated with the virtual 
motion. 

As far as forces depending on qj, gs, .-., £alone are concerned, their contribu- 
tions to P,; and Q; are evidently equal. They are different, however, in the case 
of forces depending on q,, q2, ... . Ina conservative system, for instance, the only 
forces of the last type are gyroscopic. They do no work in an actual displace- 
ment; in a virtual displacement, however, their work usually is different from 
zero. 


In contrast to its actual work, the virtual work ofa Coriolis force (Example 2. 2) 
may be different from zero. 


It follows that in a conservative system 


Q;= P;+46G;, (s=1,2,...) (5.3) 
where according to (3.1), 
_ ov oP _ oP, 
Fm bg te eee: 6.9) 


while the generalized gyroscopic forces G; are deduced from their own virtual 
work 


éW, = 2 G, 04g; (5.5) 
Thus, the alternate form of the equations of LAGRANGE, 


QoL OL : ; - 
dt Ode _ 00; = 0, (¢ — i 2, <2) (5.6) 


obtained by introduction of the kinetic potential L = T—V, is not valid, as is 
often assumed, for any conservative system. It is restricted to purely nongyro- 
scopic systems, whereas, for purely gyroscopic systems, 


ad OL OL . 
Me te (§=1.2,...) (5.7) 


Example 5.1. The particle of Example 3.1, examined in a uniformly rotating 
reference frame, has the energies 


1 : ot oor : 1 
T= zm (424+ y2+ 2%), Vu — zm wi (#8? + y%) + mee. 
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Hence, the kinetic potential takes the form 
L 1 a ies af 
ok ee oe at at: a") = 22). 


From the virtual work of the Coriolis force, 
OW, = 2m w (y dx — # Oy) , 
the generalized gyroscopic forces 
G,=2moy, Gy = —2mox, G-= 0 
are obtained. Substituting these relations into (5.7), we obtain the differential 
equations 
¥#—2oy—ow**=0, VYt+t2wo%#—o2?y=0, Ft+eg=0. (5.8) 


Using (5.6) instead of (5.7), the terms containing ¥# and 7 would have been lost. 


Under certain assumptions [12] it is possible to use (5.6) even for gyro- 
scopic systems. In these cases, however, L is no longer the difference T—V. 


6. Stability of a Configuration of Equilibrium 


In almost the whole field of physics (see for instance [17], p. 13) a configu- 
ration g, = % =... = 0 of equilibrium is called stable whenever, for any set 
Yio Vio (@ = 1, 2, ...) of initial values numerically sufficiently small, the coor- 
dinates q,(¢) and the generalized velocities ¢,(¢) remain arbitrarily small. 

In elasticity, it has become a custom to define “kinetic stability’’ in this 
way, while a system is called ‘‘statically stable” whenever the work done over 
a small displacement from the equilibrium configuration is negative. Yet there 
is no reason for using two different stability concepts. Moreover, it will be 
shown in 8 and 10 that in certain cases the two concepts are in direct opposition. 
Since, from a physical viewpoint, “kinetic stability’’ alone matters, it is reason- 
able to adopt a single stability concept as defined above. 

If 1 is the degree of freedom, configuration and state of motion of a system 
can be represented by a point (44, 9 ---» go, ..-) in the phase space of 2” 
dimensions. Thus, its equilibrium, represented by the origin O, is stable as long 
as the phase point, starting within a sufficiently small sphere of center O, re- 
mains inside a similarly located sphere of arbitrarily small radius. 

It is obvious that any stability problem can be solved by integration of the 
differential equations of motion. However, for systems containing only conser- 
vative and dissipative forces, a partial result can be established by means of 
the energy theorem. Provided that the total energy E = T+ V is a continuous 
function of q;, 9; (¢ = 1, 2,...), it can be minimized on the surface of a sphere 
g@t@g@t+...+@+@4+-::=7°* in the phase space. By physical reasons, 
T is a positive definite function of 4%, G -. Le. T=0 for g=Q=.-.=9 
and T > 0 for any other set of g,) , depending also on 4;, 9», --. - It V is a posi- 
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tive definite function of g,, g:, ..., the minimum e¢ of E on the sphere of radius 
7 is positive. Moreover, provided that the initial values q;, 9i9 are numerically 
sufficiently small, Ey<«. Since E cannot increase throughout the motion, 
E < E, < e. Hence, the phase point does not reach the surface of the sphere 7. 
We have, therefore, 

Theorem 6. If the total energy is continuous, the equilibrium of a system con- 
taining only conservative and dissipative forces is stable whenever its potential 
energy 1s positive definite. 

Theorem 6 applies to the first three categories of Table 4. It does not imply, 
however, that the system is unstable if V is not positive definite. 


II. Linear Systems 


A mechanical system shall be called /inear, if its motion is governed by a 
single set of linear differential equations. Actually, this set does not contain 
derivatives of higher than the second order; thus, it is of the form 


Dy (Min de + Bin et Cn) +4; =0, (6¢=1,2,...) (7.0) 


where the coefficients are given functions of the time. 

We shall, without mentioning it again, confine ourselves to linear, holonomic, 
scleronomic systems. Our next problem consists in the adaptation, by extending 
an idea of TH. von KARMAN and M. A. Bror [18], of (7.0) to the different 
categories of systems specified in Table 4. Hereafter it will be easy to draw 
more definite conclusions concerning the stability of every category. 


7. Purely Nongyroscopic Systems 


Purely nongyroscopic systems are entirely determined by their energies T 
and V. 


The kinetic energy of a scleronomic system is (see [19], vol. III, p. 48) a 
homogeneous quadratic function of its generalized velocities, 


1 ete 
T= > Dues (Y1 Ya +++) Gi Qe» (Mee = Mex) (7-1) 
ik 
with coefficients depending on the coordinates. For small displacements, these 


coefficients, provided they are expressible in power series, can be replaced by 
Mix(9, 9, ...) = my. Thus, T becomes a quadratic form 


1 3! 
ma 3 a Mik Vi Tk (My § Ea Mx.) (7.2) 


of the q;, defined by the constant and symmetric matrix (m,,). Since T is 
positive definite, so is (see [20], vol. I, p. 103) the matrix (m,,) (its determinant, 
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according to [21], vol. I, p. 296, being positive, together with every element of 
its chain of principal minors). 


If the potential energy V(q,, q, ...) is also expressible in a power series, 


OV 1 0?V 
WO, 0.2.) 4 re p= ——— 
oe ers 00,.Bies) Get sb 04; 99x, On DME Nae Gait tee 
(7.3) 

its constant term, which is arbitrary, can be set equal to zero. Fu#@thermore, if 
% = % = ... = 0 defines a configuration of equilibrium, 

OV : 

ae 0: = Te? Eee Fa 


hence, the linear terms, likewise, are absent. If, finally, there is at least one 
term of second degree, the coefficient of which does not vanish, V may be ap- 
proximated by the quadratic form 


Oza 


a Cin Wide |Cik = 049; O9p (0, O,...) = Crs (723) 


homed 
1,8 


ee wee 
tas 


of the q;, defined by the constant and symmetric matrix (c,,). 

Expressions (7.2) and (7.5) are approximations, restricted to terms of the 
second degree. Substituting them into (5.6), we obtain the linearized differen- 
tial equations of motion 

Dy (Mir Ge + Cex Me) =9, (@=1,2,...) (7-9) 
R 
(mx), (Czz) being constant and symmetric, (m,,) even positive definite. These 
equations represent (7.0) in the purely nongyroscopic case. 

It is clear that, in general, nongyroscopic systems are not necessarily linear. 
Nor is it always possible to linearize them. There are cases where V 1s not 
expressible in a power series; in other cases, this series contains no quadratic 
terms. 


Let m be a particle (Figure 7.1) between two springs of different stiffness ¢,, cy 
espectively. The spring force and its potential are 


i 
—(,%, |zart, (x 2 0) 
F(x) = ae 
—l,%, | Feat. (4 


IIA 


Hence, the motion is governed by either one of the two differential equations 
His oe 0 (¥% = 0) m t+ bo 4 = 0. «x S,0) 

The problem is nonlinear, since the potential V(*), the second derivative of which 

is discontinuous at 7 = 0, is not expressible in a single series. 
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If F(*) = —c x8 is the spring force of an oscillator (Figure 7.2) with a single 
degree of freedom, its potential V(%) = ¢ ¥4/4 does not contain any terms of the 
second degree. Accordingly, the differential equation of motion, 


mMm#+cxt=0, 


when linearized, does not represent an oscillation. 


Cp Cy F 
x 
m : ‘ 
F x 
x 
Figure 7.1 Figure 7.2 
Nonlinear oscillator. Nonlinear oscillator. 


Since (7.2) is positive definite, there exists a linear transformation 
Pim Dd) Gin Me (¢=1,2,...) (7.7) 
k 
(see [22], vol. I, p. 32, or [23], p. 188), which yields purely quadratic forms 
lg 1 °°. y 1 9 
Sabre Ci a Tt (7.8) 


for both energies. The quantities q,, @,..., vanishing in the equilibrium con- 
figuration, are the normal coordinates of the system; (7.8) are the normal forms 
of T and V. The coefficients m,, c; are real; the m,’s are even positive. Substi- 
tuting (7.8) into (5.6), we obtain the differential equations 


mGitceg=0, ($001, 2,405) (7,9) 


each of which contains no more than one normal coordinate. 
The differential equations (7.9) are solved by setting 


g=A,e%, B= —4, (¢=1,2,...) (7.10) 


the A; being arbitrary constants. In general, (7.10) represents two fundamental 
solutions for every normal coordinate, characterized by the roots A;, —A;. Pro- 
vided the coefficients ¢; are positive, all exponents are imaginary (Figure 7.3). 
In this case, the fundamental solutions are harmonic oscillations. If at least 
one ¢; is negative, the corresponding roots are real (Figure 7.4). Since one 
of them is positive, there exists an illimited fundamental solution. If c, = 0, 
(7.10) yields but one fundamental solution »; = A, for the 7 in question, the 
other one being p; = B;, t. It follows that the system is stable only as long as 
all coefficients c; are positive. Since c; > 0 (i = 1, 2, ...) also are the conditions 
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necessary and sufficient in order that the matrix (C:,) Or, equivalently, the 
quadratic form V be positive definite, Theorem 6, for the case in question, can 
be restated as follows: 

Theorem 7. The equilibrium of a purely nongyroscopic system ts stable only so 
long as its potential energy is positive definite. 


Figure 7.3 Figure 7.4 Figure 7.5 
Characteristic roots of a purely Characteristic roots of a purely Particle tied elastically to the 
nongyroscopic system nongyroscopic system origin. 
(V being positive definite). (V having negative values). 


Example 7.1. The differential equations governing the motion of a particle m 
(Figure 7.5) tied elastically to the origin are 


Mae Cr=O0, mMy+cy=0, mztcz=0. (7.11) 


Here, %,y,z already are normal coordinates. The equilibrium position 
y= y = = 0 is stable, provided c>0. 


8. Purely Gyroscopic Systems 


According to 2, both gyroscopic and dissipative forces depend on the genera- 
lized velocities. Thus, purely gyroscopic and purely dissipative systems are 
characterized by generalized forces of the form 


G,= —D)) gin de» (scent, 26 lel S) 
k 


the coefficients g;; being given functions of ¢. Their differential equations of 
motion, instead of (7.6), are 


Dy (Mine t+ ix Ie + Cin Me) = 9- ae oa 7a a (S794) 
D 


The power of the forces (8.1) is 
P=) G,q¢=—D) in 1% - (8.3) 
4 1,k 


If (gix), (g/,) denote the symmetric and the antisymmetric parts of the 


matrix (g;x), _ 
P=—)) Bindi de (8.4) 
i,k 


104 Hans ZIEGLER ZAMP 


since the sum containing (g/,,) is zero. According to 2, gyroscopic forces do no 
work in actual motions, whereas the work of dissipative forces is negative. 
Hence, the matrix (g/,,) represents gyroscopic forces!), while (g;,), provided it 
is positive definite, represents dissipative forces. Thus, the differential equations 
of motion of a purely gyroscopic system, in particular, are given by (8.2), 
where (m;,) and (c;,,) are constant and symmetric. (m,,) is positive definite and 
(g;x) 1s antisymmetric. 

According to Theorem 6, the equilibrium of a purely gyroscopic system is 
stable as long as V is positive definite. By comparison with Theorem 7 follows 

Theorem 8a. A conservative system cannot be made unstable by gyroscopic 
forces. 

In the case of a constant matrix (g;,) the differential equations (8.2) are 
solved by setting 


gi = A,e™. (4.2 1,2, 25.) 0 ofan 
Substitution of (8.5) into (8.2) yields the characteristic system 


Dy (mind? + in A+ Cix) Ap =O. (¢=1,2,...) (8.6) 
k 


Excluding the trivial solution A, = A, = ... = 0, we obtain the characteristic 
equation 


|| mp A? + Bin A+ Cy || = 0. (8.7) 

The value of the determinant in (8. 7) is not affected if rows and columns 
are interchanged. Due to the symmetry of (m,,), (c;,) and to the antisymmetry 
of (g;,), such an interchange is equivalent to a change of sign in 4. Hence, if A; 
is a solution of (8.7), so is —A,. Since a positive real part in any root means 
instability, the system is stable only as long as all roots are imaginary 
(Figure 7.3). This condition, according to Theorem 6, is satisfied when V is 
positive definite. . 


Thus, Theorem 6 implies that, if (mx), (€;,) are Constant, symmetric and posi- 
tive definite, and if (g;,) is constant and antisymmetric, (8.7) has only imaginary 
roots. 


If, on the other hand, V is not positive definite, the system, according to 
Theorem 6, would be unstable in the absence of gyroscopic forces (Figure 7.4). 
There is, however, no reason why it should be unstable in their presence. 


Example 8.1. Let m be a particle (Figure 8.1) tied elastically to the axes of 
a uniformly rotating plane coordinate frame. Its energies are given by 


et : eet 
T= =m (#7? +4"), J = 3 (a — mw?) H24 (Cg — mo?*) y*), 


4 


1) See also [40], vol. I, p. 258. 
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the gyroscopic forces by 


Gz=2moy, Gy=—2mo#; 


hence, the differential equations of motion are 


Be é 
#-209+(2 wt) =0, YVt2wke ay = 
Yt2Zuoe+ shake 5) == (0). (8.8) 


WwW 
y mary 
CX mw?x 
m > ‘ 
io) 

oy A 
(a) SS 
2M Wx x 

Figure 8.1 


Relative motion under elastic forces. 


The characteristic equation is 


C Cc G C 
Af ee eed ie, eee 202 py as heh 2 2 pe Q 
+ (a+ + wr) A le wt) (58 o) (ON (8. 9) 
Obviously, its roots 
1 C 6 1 = 2 
hie om a : oF |/3 (a) eee oe 


provided c,>c,>0, are real and of different absolute values. In consequence, 
there are four roots 4, real or imaginary, all different if none of them are zero. 
Thus, if c,/m+ w?+ c,/m, (8.5) yields four fundamental solutions, and the 
general solution is limited so long as all values of A are imaginary, i. e. so long as 
both 4? are negative. 

If w2 <c,/m <c,/m, the coefficients in (8.9) are positive. Therefore, 1? < 0, 
73 < 0: the system, in accordance with Theorem 6, is stable. For ¢,/m < w? <c,/m 
the potential energy V is no longer positive definite. The constant term in (8.9) 
now is negative; thus, Aj < 0, 4} > 0: the system is unstable. If ¢,/m <c./m < w?, 
both roots 22, 42 are again negative: the system is stable while V is negative 


definite. 


Under the last of these assumptions Figure 8.1 represents @ system which, 
according to 6, is ‘‘kinetically stable” though “statically unstable’. 

From this example follows 

Theorem 8b. If the potential energy of a purely gyroscopic system 1s not 
positive definite, the equilibrium may be stable or unstable. 

By comparison with Theorem 6 follows 

Theorem 8c. A conservative system may be stabilized by gyroscopic forces. 
ses w2 =¢,/m <c,/m and 
m being tied elastically 
>c/m meet. Since for 


It is easy to see that in Example 8.1 the limiting ca 
c,/m <¢,/m =o? are unstable. If, however, cy = =¢ ( 
to the origin), the stable domains w2<c/m and w? 
w2 = c/m, the differential equations are 


59 a) ay Gh ey. y+t2Zo#=0, (8.10) 
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having the general solutions 
*—Acos2wit+ Bsin2Zwit+C, y=-—A sin2wt+Bcos2wt+D, 


the equilibrium position 7 = y = 0 is always stable. For an observer taking part 
in the rotation, the system appears to be stabilized for m* = c/m by the Coriolis 
force. For an observer at rest, the result of course is evident. 


9. Purely Dissipative Systems 


According to Theorem 6, the equilibrium of a purely dissipative system is 
stable as long as V is positive definite. By comparison with Theorem 7 follows 
Theorem 9a. A purely nongyroscopic system cannot be made unstable by — 
dissipative forces. 


According to 8, the differential equations of motion of a “, 
purely dissipative system are (8.2), where (mx), (Six), (Cix) Yi 
are symmetric, (m;x), (cix) are constant and (mx), (gx) are 
positive definite. 

In a conservative system the gyroscopic forces may be 
functions of ¢t. However, dissipative forces depending on ¢, 
according to Table 4, would render the system instationary. 
Thus, in a purely dissipative system (g;,) is constant. Hence, ph 
the differential equations (8.2) are solved by (8.5), the char- Chane 

Te : : ; roots of a purely 

acteristic equation still being (8.7). Its roots A are real or dissipative system 
conjugate complex. Conjugate imaginary roots (as in Figure _(V being positive 
7.3), representing undamped oscillations, appear no more. definitely 


According to Theorem 6, none of the roots have a positive 
real part while V is positive definite. Thus, every root is shifted, by dissipation, 
from the imaginary axis to the stable half-plane (Figure 9.1). 


Theorem 6 thus implies that, provided (m,,), (g;,), (¢;,) are constant, sym- 
metric and positive definite, all of the roots of (8.7) have negative real parts. 

Example 9.1, If the particle of Example 7.1 is subjected to a damping force 
R= --2myv, its differential equations of motion are 


MB 2: HANG, HOV MN hee (9.1) 
The characteristic equation, 


A+ 2y Ae SO, 
m 


A= +» #|/n-2 i . 
We 


Tor clm> 0 both roots have negative real parts: the system, in accordance with 
Theorem 6, is stable. If c/m <0, one of the roots is positive: the system, as in 
the absence of damping, is unstable. For c/m = 0 the differential equations 
reduce to 


has two roots 


H+2yeH0,..5 00. (9. 2) 
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Their general solution is 
= Bin, 
MA Wo gO. OP cr inn 


Thus, by contrast to the undamped case, the equilibrium position ¥ = y = z= 0 
is still stable although V has ceased to be positive definite. 


It follows from this example that Theorem 8b is also valid for purely dissi- 
pative systems. To be sure, the damping force in Example 9.1 stabilizes the 
system only in the limiting case when V is just losing its positive definite 
character. It is, however, more important to know the effect.of small dissipa- 
tive forces beyond this limit than to generalize the last result. 

Provided that one of the coefficients c; in the normal form (7.8) of V is 
positive, the positive real axis in Figure 7.4 contains, in the absence of dissi- 
pative forces, at least one root. Since it is a continuous function of the coeffi- 
cients g;,, it remains outside the stable domain while the dissipative forces are 
sufficiently small. We have, therefore, 

Theorem 9b. If, in a purely dissipative system, there exists a configuration in 
which V <0, the equilibrium is unstable, provided the dissipative forces are suffi- 
ciently small. 

By comparison with Theorem 7 we obtain 

Theorem 9c. Whenever a configuration exists in which V <0, a purely non- 
gyroscopic system will not be stabilized by sufficiently small dissipative forces. 


10. Purely Circulatory Systems 


Since, according to 2, circulatory forces depend on the coordinates only, 
purely circulatory systems are charcterized by generalized forces of the form 


fee = Ode, ye ite han aren 8 
k 


the coefficients c,,, being constant. 
If the matrix (c,,) were symmetric, the forces (10.1) might be derived from 


a single-valued potential 
1 
Vez Did Tk + (10.2) 
Pe Gok 3 


In consequence, the system would be conservative. Thus, (¢;,) is asymmetric ; 
its symmetric part (c;,) represents the conservative forces already accounted 
for in 7, while its antisymmetric part (cj,) represents the circulatory forces. 
Again, the differential equations of a purely circulatory system are (7.6), 
where (m;,), (Cix) are constant, (m;,) is symmetric and positive definite, and 
(c;,) is asymmetric. They are solved using (8.5). The characteristic equation 1s 


|| men A® + Cex | =9. (10.3) 


If 2, is a solution of (10.3), so is —j,. Hence, the equilibrium is stable exactly 
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so long as all roots A; are imaginary (Figure 7.3). In the absence of circulato 
forces this condition, in accordance with Theorem 7, is satisfied only so long as 
V is positive definite. Yet, if circulatory forces are present, there is no reason 
for the system to behave in the same manner. 

Example 10.7. The particle m of Figure 10.1, tied elastically to the axes 
a fixed plane coordinate frame, and subjected to a force of magnitude yw /, actin 


perpendicular to the radius vector, is circulatory. Its energi 
are given by 


oye 
Ta=>m(#+9), Vazlaxtt+ ay), 


and the generalized circulatory forces by 


F,=—py, R=pe. 


Figure 10.1 
Elasticforcesand force Lhe differential equations of motion are 


perpendicular to the 
Avaric es c L 
radius. pina eget ae y=0 y+ Bet A ale, 4 (10. 


They yield the characteristic equation 


2 
q+ Ce 42 , 16 + 


24 \ _ 0) 
ee m § m2 
the roots of which are 
1 = 
Je eth ge" (cy + eg + V (eq — &)*? — 4 w*) 
lf ¢,=@=c>0, 
1 
Aa = — — L 

ule Et B): 


thus, the system is unstable while its potential energy is positive definite. 
= 2, C= —p/3, 
me ee) oe 


6 m ’ 


hence, the system is stable while its potential energy is not definite. 


Under the first assumption, Figure 10.1 represents a system which, accor 
ing to 6, is “kinetically unstable” though “statically stable’. Under the ad 
assumption, the reverse is true. 
From this example follows 
Theorem 10a. The equilibrium of a purely circulatory system can be sab 
or unstable regardless of tts potential energy being positive definite or not. 
By comparison with Theorem 7 we obtain 
Theorem 10b. Conservative systems may be stabilized or made unstable / 


ee 


circulatory forces. 
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11. Constant Generalized Forces 


The properties of the systems treated thus far do not alter if constant 
generalized forces —h,, —hy,... are present. Their differential equations, 
however, become inhomogeneous. 

The general solution of (7.0) is obtained by combination of the general 
solution of the reduced system (8.2) with a particular integral of (7.0). The 
first has been treated in 7 to 10; it represents the motion of the system in the 
absence of constant forces. The second may be taken as 


G=A;, (grata Deeie NS (11d) 


where the constants A, are the solutions of 


D>, Cx Ar == —h,. (2 = Il, Zs eA (alla 

This represents the equilibrium configuration which, under the action of con- 
stant forces, is no longer q, = g, =... = 0. If new coordinates 

Dec ean (@=1,2,...) (11.3) 


are introduced into (7.0), we obtain, in view of (11.2), the homogeneous dif- 
ferential equations 


Dj (Mex Put Bin be + Cox Pe) = 9. (Soe bora roliblecs) 


They are equivalent to (8.2) and represent the motion of the system with 
reference to the equilibrium configuration (11.1). 

The system is endangered whenever the general solution of the reduced set 
(8.2) is illimited, i.e. whenever the equilibrium (11.1) is unstable. Thus, the 
cases of instability examined in 6 to 10 are dangerous, whether constant forces 
are present or not. 

Another danger, however, arises whenever the particular integral (11.1) 
becomes illimited, i. e. whenever the equilibrium configuration, under the ac- 
tion of the constant forces, becomes infinitely removed. This phenomenon, 
which is independent of the stability of (11.1), would pass unnoticed if the 
constant forces were disregarded; it may, therefore, be called a latent insta- 
bility of the system. 

A latent instability arises whenever the determinant | cz, || of (11.2) 
vanishes. According to (8.2), however, in the absence of constant forces, 
|| c;, || = also is the condition necessary and sufficient for the existence of 
(an infinity of) equilibrium configurations apart front G) = ¢ =. = OF Lm 
order to distinguish these additional configurations of equilibrium from those 
considered so far, we call them nontrivial equilibrium configurations, we have 


; 
, 
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Theorem 11a. In conservative, purely dissipative and purely circulatory sys- 
tems the stability of the equilibrium remains unaffected by constant forces. There 
appear, however, latent instabilities whenever the system, in the absence of constant 
forces, has nontrivial equilibrium positions. 

Evidently, Theorem 11a remains valid if conservative, dissipative and circu- 


latory forces are present at the same time. It is even valid if the matrices (m,,), 
(g;;,) depend on the time, i. e. for special classes of instationary systems. 


Example 71.7. If, in Example 10.1, in addition to the forces shown in 
Figure 10.1, the particle is subjected to constant forces X, Y (Figure 11.1), the 
differential equations of motion become 


ek el oe eee n ¥, . 
Pe ye ey eee. fh tN gE ag eee (11.5) 
m m - m m m m 4 
The trivial equilibrium position is 
. CoX — eh Y el Seat or 45, 
SS byte ee Cy 6 + pe? 


If cy cg + w? = 0, it is infinitely removed; at the same time, in the absence of AG 
Y, there appear nontrivial equilibrium positions 


S48 acme ae } 

C - le Cy ; 

: 

in addition to the trivial one, x = Vyie=wOe : 


In the absence of circulatory forces the determinant | Con || is symmetric, 
If normal coordinates (7.7) are used, || c;, || = cy Cy cs .... Thus, we obtain : 

Theorem 11b. In conservative and purely dissipative systems there exist non- 
trivial equilibrium configurations whenever at least one coefficient c; in the normal 
form of V vanishes. 


Obviously, Theorem 11b remains valid if conservative, dissipative and even 
instationary forces are present simultaneously, provided the last ones appear in 
the matrices (m,,), (g;,) only. ; 


Suppose the coefficients c;, originally all positive, decrease. As soon as one 
of them becomes zero, the first nontrivial equilibrium configurations appear. 
Simultaneously, V ceases to be positive definite. Thus, the limit which, accord- 
ing to Theorems 6, 7, and 9b plays the decisive role in the stability of conser- 
vative and purely dissipative systems, is marked by the first appearance of non- 
trivial equilibrium configurations. 


In Examples 7.1 and 9.1, for c= 0 any point in space represents a position 
of equilibrium. 
Example 77.2. Vf, in Example 8.1, in addition to the forces shown in Figure 8.1, 


constant forces X, Y are applied (Figure 11.2), the differential equations of 
motion are 


# fC xX + c yA 
pra D 1 2 Syhakt y Y 2 2 — 
x ay +| wt) x ae | ow) y - (11.6% 


Vol. IV, 1953 Linear Elastic Stability itil 


Provided c, < Cy, the trivial equilibrium position 


xX Wy 
4 = ——_—,, vy = ———— 
Cy — Mw? Co — mw? 
is unstable while ¢,/m < w? < ¢,/m. The limits of this interval mark the existence 
of nontrivial equilibrium positions (in the absence of X, Y) and, at the same 
time, of latent instabilities (in their presence). 


1a 
¥ meviy | 
moe x Xx 
—n 
Gx A 2mey 
Cau, 
@ )- [ > 
emax 
Figure 11.1 Figure 11.2 
Particle of Figure 10.1 acted upon by Particle of Figure 8.1 acted upon by 
constant forces. constant forces. 


If cy = Cy, the equilibrium, due to the stabilizing effect of the Coriolis force, 
is stable for any value of w; however, the latent instability at w? = c/m is still 
present. 


12. Purely Instationary Systems 


If in (7.0) the coefficients depend on the time, the system, provided that 
not only the gyroscopic terms are functions of #, is instationary. In order that 
it be purely instationary in the sense of Table 4, (7.0) must have the form 


DS (min tet ng) += 0, (=1,2,...) (21) 
& 


(m;,,) being constant, symmetric and positive definite while (c;;) 
and h; are given functions of ?. 
It is clear that systems of this kind do not necessarily have 


an equilibrium configuration. If one exists, there is no reason 


why its stability should obey any of the laws established so far. le 
Mp Cos 
Example 12.1. Let m (Figure 12.1) be a mathematical pendu- gat 
lum, subjected to its weight anda harmonically oscillating pertur- Figure 12.1 
bation. The differential equation of motion, when linearized, is Pendulum with 
} oscillating 
G+ (4+ Foose i) = 0. (12.2) axial load. 


If we let t = wt, a new time scale is introduced. By designating derivatives with 
respect to t by primes, we have 

6” + (6 + ecost) B= 0, 
where 
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This differential equation, which is of the MATHIEU type, may be deduced from 
the energies 


(0 -+ ecost) 8; 


since V depends on t, however, the problem is nonconservative. 
The general solution is obtained by superposition of two fundamental solutions ; 

whether it is limited or not depends on the parameters 0, ¢, i. e. on the period and 

the intensity of the perturbation. Figure 12.2 (taken from [24] and often referred 


i\ 1 JAC 
may 

| eC . 
PANINNSS KN SG 


Figure 12.2 
Strrutt’s Table. 


to as Strutt’s Table) illustrates the stable and unstable sets 6, e. The shaded 
regions, which are symmetric with respect to the 6-axis, represent the sets for 
which the motion is stable. Any point outside the shaded regions (boundary 
points with the exception of those on the d-axis included) gives rise to an unstable 
solution. 

In the absence of a perturbation, the representative point P(6, e) in Figure 12.2 
lies on the positive 6-axis; thus, the pendulum is stable. When # increases, P moves 
parallel to the e-axis. No matter what the value of 6 is, P can, by a suitable 
choice of #, be made to leave the shaded region. Thus, it is always possible to make 
the pendulum unstable by a perturbation of any given period. 

If / is assumed to be negative, the inverted pendulum is obtained. Its repre- 
sentative point P, in the absence of a perturbation, lies on the negative d-axis; 
hence, the pendulum is unstable. However, for any value of 6 by a suitable 
choice of ~ the point P can be shifted into the interior of a branch of the shaded 
region; thus, the pendulum becomes stable. 

By the straight line e = 6 the upper half plane of Figure 12.2 is divided into two 
regions I, II. When P lies inside I, V is positive definite and vice versa. 


From this example follows 

Theorem 12a. The equilibrium of a purely instationary system can be stable 
or unstable regardless of its potential energy being positive definite or not. 

By comparison with Theorem 7 we obtain 


Theorem 12b. Conservative systems may be stabilized or made unstable by 
instationary forces. 
‘ 


a a iene 
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III. Elastic Stability 


By means of the classification of Table 4 and the theorems established in 
6 to 12, we are ready to analyze the different methods in use to solve stability 
problems, particularly in regard to their limits of validity. 


13. Methods for the Calculation of Critical Loads 


In the unloaded state any elastic system has a configuration of equilibrium 
which is stable provided the system is not prestressed. If under the action of 
a given load and of small though arbitrary initial disturbances, the displace- 

ments qj, 9, -.. measured from this equilibrium position become too great, the 
system is endangered. 

Ordinary problems of linear elasticity are characterized by small displace- 
ments proportional to the applied load. If these displacements, under certain 
values of the load, increase disproportionately and without a limit, we encoun- 
ter a stability problem. Here the increase of the quantities q,, qs, ... may be 
due either to the displacement of the equilibrium configuration itself or to 
its destabilization by the load. 

It is obvious (see 6) that any stability problem can be solved by integration 
of the differential equations of motion, i.e. by the kinetic method. The loads 
for which, at least under the influence of small initial disturbances, illimited 
solutions exist, are called critical loads. Thus, we have 

Theorem 13: The kinetic method is valid without restriction. 


Evidently, it even applies in nonlinear, nonholonomic and rheonomic cases. 


In rotating systems, at least some of the loads are centrifugal forces. In 
such cases, every critical load is associated with a certain angular velocity; 
thus, we are faced with a problem of critical speeds. 

Usually the application of the kinetic method is cumbersome. Therefore, it 
is not surprising that it has been avoided as far as possible, and that, for sta- 
tionary systems, the simpler static methods are preferred. 

One of them is the equilibrium method already used by L. EULER in his 
famous investigation of compressed rods [25]. It defines as critical the values 
of the load for which nontrivial equilibrium positions exist. 

The other one is the energy method developed by S. TIMOSHENKO (see [1], 
p. 78) and others. It defines as critical the value of the load for which the 
potential energy of the system ceases to be positive definite. 

In a slightly different interpretation, the energy method defines as critical 
the values of the load for which the potential energy is stationary; in this modi- 
fied form the method is equivalent to the equilibrium method. 

Both static methods have been applied with such succes that neither their 
validity nor their equivalence with the kinetic method ever have been questioned. 


ZAMP IV/8 
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They are used so exclusively that the concept of stability (which, according to 6, 
is essentially kinetic) in elasticity has acquired a purely static meaning (which, 
according to 8 and 10, occasionally results in errors). 


: 
14. Purely Nongyroscopic Systems 


Most buckling problems are of the purely nongyroscopic type. The poten- 
tial energy of the system is the sum 


V=U+W (14.1) 
of the strain energy U and the potential energy 
W=-PW, (14.2) 


of the load P [8]. Whereas U is positive definite, W, is not negative definite; 
hence, the coefficients c; in the normal form of V, positive for sufficiently 
small values of P, decrease with increasing P and (at least part of them) 
change their sign. Under the load P, for which the first of them, c,, passes 
through zero, V ceases to be positive definite; at the same time, by Theorem 7, 
the equilibrium becomes unstable (Table 14). Moreover, whenever c; = 0, i. e. 
under a series of loads P, < P,<..., there appear, according to Theorem 11b, 
nontrivial equilibrium configurations. Thus, we obtain at once 


Table 14 
Critical loads for systems of different categories 


V not pos. definite 


Zi 


nontrivial equilibrium positions 


V pos definite 


2 
Pf fe 
Vy, ie 


stable regions unstable regions stability depending 
on case 


Vol. IV, 1953 Linear Elastic Stability MUS 


Theorem 14a. In purely nongyroscopic systems, the energy method yields the 
smallest critical load P,. The equilibrium method, too, yields P, as the smallest 
load of a series P< P,<-:-, the other elements of which are insignificant, 
since any load P = P, ts critical. 


Example 14.1. Figure 14.1 shows a simplified model of an elastic rod built 
in at one end and compressed by an axial load. The potential energy, 


is positive definite as long as P< c/l. The equilibrium condition is 
G—PiIb=0, 
and the differential equation of motion is 


o + ee P20 (14.3) 


Thus, by either of the static methods as well as by the kinetic method the critical 
loads P= P, = c/l are obtained. 


Actual systems never are free from physical imperfections. As a consequence, 
in stationary systems, the trivial equilibrium configuration is slightly different 
from 9,=4.=.-.=0. Any such deviation, according to 11, is equivalent to 
the presence of constant forces. Therefore, the possibility of latent instabilities 
must be considered. However, according to Theorems 11, latent instabilities 
appear whenever c; = 0, i.e. under the loads P,, P,... for which the equilib- 
rium itself is unstable. Thus, we have 

Theorem 14b: In purely nongyroscopic systems, latent instabilities due to 
physical imperfections can be entirely disregarded. 


8 3 
cd , cob z 
o\p 


7) 
M 
Figure 14.1 Figure 14.2 
Pendulum loaded axially. Pendulum of Figure 14.1 acted upon by constant forces. 


Example 14.2. lf, in Example 14.1, the load is slightly oblique and eccentric, 
these imperfections can be represented by the small constant quantities Q, M 
(Figure 14.2). The equilibrium condition is 


(c-P)6-QI-M=0. 


-e. when the smallest 


Thus, the only latent instability arises when P= he aera 
critical load of the perfect pendulum (Q = M = 0) is reached. 
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15. Purely Gyroscopic Systems 


Most problems of critical speeds (if referred to a rotating coordinate frame) 
are of the purely gyroscopic type. According to Theorems 8, they differ from 
those considered in 14 insofar as the gyroscopic forces are capable of stabilizing 
the equilibrium for P = P,. By Theorems 11, however, they do not stabilize 
the latent instabilities under the loads P,, P,,... (Table 14). Thus, we have 

Theorem 15a: In purely gyroscopic systems, the energy method yields the 
smallest critical load P,. The equilibrium method yields the critical loads 
P,<B<....Any other load P> P, may or may not be critical. 

Thus, by the static methods alone it is impossible to decide whether a load 
P> P, but not equal to B,, P,, ... is critical or not. 

In Example 8.1, provided c,>c,, any angular velocity satisfying the ine- 
quality c,/l < w? <c,/I is critical. However, the energy method only yields 
? = ¢,/l, i.e. the smallest value, whereas by means of the equilibrium method 
both limits mj = c¢,/l, w3 = c,/l of the critical interval are obtained. That the 
whole interval is critical follows only from the analysis based on the kinetic 
method. 

By contrast with 14, the loads P,, BR, ..., due to the stabilizing effect of 
the gyroscopic forces, may represent but latent instabilities which might 
remain undetected if physical imperfections were disregarded. We have, there- 
fore, 

_ Theorem 15b. In any but purely nongyroscopic systems, the existence of 
latent instabilities due to physical imperfections must be recognized. 

If, in Example 8.1, cy = cy, the equilibrium itself is stable for any angular 
velocity. For w? = c/l, however, as shown in Example 11.2, a latent instability 
arises. 

Examples 8.1 and 11.2 represent a simple problem of critical speed. For 
a shaft carrying a single disc the forces in Figure 8.1 exemplify the elastic 
restoring force, the centrifugal and the Coriolis forces for an observer taking 
part in the rotation. The additional forces X, Y in Figure 11.2 represent a 
small eccentricity of the disc. If the shaft has two equal flexural rigidities, 
Cy = Cy). Then the critical angular velocity @, is due to a latent instability and 
would escape detection if the eccentricity were neglected. 


16. The Influence of Damping 


Mechanical systems always are subjected to damping forces which, however, 
usually are small. From Theorems 9 follows (Table 14) 

Theorem 16a. Small dissipative forces do not affect the stability of a purely 
nongyroscopic system, 


There is, however, no reason to assume that other systems, likewise, are 
insensitive to damping. 
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Example ON ( ey Lik pein Example 8.1, in addition to the forces shown in Figure 8.1, 
a damping force —2 mb v is introduced (Figure 16.1), the differential equations 
of motion are, provided c, = c, = c, 
ene : 5 C 
H+260% —2m0y + (j-m a= 0, V+2Zoxv+2bdy4+ (= 08) y = Om 
m ‘ 
(16.1) 


Figure 16.1 
Particle of Figure 8.1 subjected to damping. 


The characteristic equation is 


At pi A pz A” bPsA+ fy= 9, 


where 


c Cee: he AN? 
p=4db, p= 2 (2 +o? + 204), py= 46(—— 0%), p= (Zo). 


In order that none of the four roots have a positive real part, the coefficients of 
the characteristic equation must satisfy the conditions of E. J. RourH (see [15], 
wollen tiesp. 228), 


(Oy = 0 = 0), 
pr Pr — P= 40( + 3.08 4 408) Si, 


(1 Pa — Ps) Po — Pi Dy = 64 (5 — 0?) (w* + D4) > 0, 
p= (5-07) Soe 


Evidently, they are satisfied exactly as long as w? < c/m = w{. Thus, for o>a,, 
i.e. in the range where the particle, according to Example 8.1, is stabilized by 
the Coriolis force, it is again made unstable by the damping force. 


From this example follows 

Theorem 16b. Dissipative forces, applied to other than purely nongyroscopic 
systems, may have a destabilizing effect. 

In accordance with the statement made at the end of 15, Example 16.1 
represents the problem of critical speed for a shaft equipped with a single disc 
and subjected to internal damping. The result obtained implies that any angu- 
lar velocity w >, should be critical. It contradicts the fact that, actually, 
cw, alone is critical while, for w > @y, shafts of this type run smoothly. Tt is 
improbable that this discrepancy is due to the fact that external damping 
(air resistance) has been neglected. Equally unsatisfactory results obtained in 
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other cases (see [8], p. 55) suggest that, at least in other than purely nongyro- 
scopic problems, internal damping is inadequately represented by linear damping 
forces. As a matter of fact, internal damping is characterized by a time effect 
(see [26], p. 237) which is suppressed if the damping forces are assumed to be 
linear combinations of the generalized velocities. Thus, it is to be expected that 
the actual behaviour of other systems, too, is more satisfactorily explained by 
neglecting damping: 

Theorem 16c. While in purely nongyroscopic systems it is allowed to disregard 
small damping forces, in any other system it is to be recommended. 


17. Circulatory Systems 


Certain elastic stability problems are of the circulatory type. In view of 
Theorems 10 it is not to be expected that they can be solved by static methods. 


Example 17.1. The double pendulum of Figure 17.1, characterized by concen- 
trated masses and elastic restoring moments, is a simplified model of an elastic 


Figure 17.1 


Double pendulum, loaded tangentially. 


rod built in at one end and compressed by a tangential end load. Its energies are 
i i ee sos : el 2 1 
T= zim ai 8} + m, (13, + a, %)2], V= => ¢ [8] + (8, — 4,)*), 

the generalized forces are 


Pw Pid, -0), P= 0. 


Thus, we obtain the differential equations of motion 


(m, a? + m, 1?) 6, + ml a, 6, +(2c— P )da—(c-—P1IOd=0, 
my 1 dy O, + my ad B, —c#, +cB,=0. 
The characteristic equation is 


po At + P,A*+ p= 0, 


(17.1) 


where 
Po = MyM, a5 as, 
P, = (m, ai + m,1*) c+ m, a, (1+ ay) (2c— Pl), 


P,=c*. 
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Since p, > 0, there does not exist any nontrivial equilibrium position. By Theo- 
rem 11a, latent instabilities are absent. Moreover, neither the equilibrium method 
nor the energy method yield a critical load. 

Let m,= m,=m and a,= a,=1/2. Then the masses are concentrated at 
the two centers of gravity, and Figure 17.1 represents a prismatic, homogeneous 
elastic rod. The coefficients of the characteristic equation become 


1 1 
Pag Pa eee SILO) ee Goi ae 
and its discriminant is given by 
3 
A= p} — HY alae ads eae (30-226 FP b+ 3 P2]3)" 


Provided P is sufficiently small, the quantities A, po, py, py are positive; the roots 
2? are negative and thus the pendulum is stable. When 


A is negative; the roots A? are conjugate complex: the system is unstable. If P 
exceeds the upper limit, A is positive again; however, since p, is negative, both 
roots A? are positive: the system remains unstable. Hence, any load P> 3 c/l is 
critical. 

Example 17.1 illustrates the possibility of a complete failure of the static 
methods when applied to circulatory systems. 

Wetin Example 17.1, m,= 2m, m,—m, a,=a,—=/. Thus, the mass is conc¢en- 
trated in three points. Now we obtain 


Dig Lo PAT OR ie Pig Cor, 
and 
A= m?l4 (41 c? — 28c P1l+ 4 Pl’). 


Again, the pendulum is stable for sufficiently small values of P. When 
7 —\ 6 7 NG 
(e> y2) at = fee V2) > 
A is negative, and for higher values of P the coefficient p, becomes negative. 


Thus, in contrast to the former result, any load P> (7/2 — y2) c/l = 2-09 ¢/f is 
critical. 


The instabilities appearing in Example 17.1 have the character of self-excited 
oscillations, depending on the mass distribution. Since, in static investigations, 


the masses do not enter, we have (Table 14) , 
Theorem 17: Stability problems of the circulatory type cannot be solved by 


static methods. 
18. Purely Instationary Systems 


In view of the results obtained in 17, it would be surprising if stability 
problems of the instationary type could be treated by static methods. 


: 
; 
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Example 18.1. The pendulum of Figure 18.1 represents a simplified model 
of an elastic rod built in at one end and compressed by a pulsating axial load. 
Its differential equation of motion is 


ne futoeat ta. eéeds ‘) o— 0. (18.1) 
Taking, asin Example 12.1, r= /, and introducing the notations 
Rp ei ois ee ee Pe hii 
m a®w?’ m a? wm?’ 
Figure 18.1 we obtain the Mathieu equation 


Pendulum 
under pulsating 
axial load. 


0” + (6 +ecost) B= 0, (18. 2) 


the solution of which is stable whenever in StRut?T’s Table (Fig- 
ure 12.2) the point P(6, «) lies in the shaded region. 

If the mass of the pendulum is concentrated at its center, a = //2, 

4c d 4p 


, &= : 
m 1? w? ml ow? 


If, on the other hand, it is distributed between the two end sections, we obtain, 
replacing m by m/2 and taking a=1, 
Ac 1 2 


She ee ae : eal iS ot 
m 12 w? 2 milo? 


As is seen in Figure 12.2, it is possible that of the two points P(6d, e), P’(d/2, ¢/2) 
only one is stable. Thus, here too, stability depends on the mass distribution. 


From this example follows (Table 14) 


Theorem 18. Stability problems of the purely instationary type cannot be 
solved by static methods. 


This theorem, of course, does not apply to the more general case where gyro- 
scopic or dissipative forces alone are instationary. 
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Untersuchung einiger Integrale mit Bessel-Funktionen, 
die fiir die Elastizititstheorie von Bedeutung sind 


Von JOHANNES DOrR, Darmstadt?) 


1. Aufgabenstellung 


In neueren Untersuchungen der Elastizitatstheorie®) treten Integrale der 
Form 


rs ym 2(f¢ 
Ti) (1, t,) -/% z ae ay dt, (m=0,1,2) 
0 


auf, worin die ¢, und ¢, Nullstellen der Bessel-Funktion /,(¢) sind. Durch 
Partialbruchzerlegungen lassen sich diese Integrale auf die folgenden vier 
Typen zuriickfihren: 


Psi) 20 ee Tid) 
Ot t,) —— ioe th at , R(+ 9) > (t a3 4,3 dt . 
0 0 
Die Integrale R(-{ ¢,) treten nur im Falle k = » auf. Die numerische Berech- 
nung dieser Integrale ist wegen des unendlichen Integrationsweges miihsam 
und ungenau. Ein eleganterer, von SZABO angegebener Weg bietet sich durch 


Reihenentwicklung des Integranden, indem man ¢ + ?¢, = z setzt und auf die 
Reihe 


(2 t,) = cs - Ii Ai 2) Jin( i t,) 


zuriickgreift. Durch Quadrieren erhalt man daraus eine Doppelreihe, mit deren 
Hilfe sich die Integrale Q und R berechnen lassen, wobei allerdings der Arbeits- 
aufwand wegen der Doppelsummen recht betriachtlich ist. 

Im folgenden wird ein Weg gezeigt, die Integrale mit Hilfe numerischer 
Integration (zum Beispiel mit der Simpson-Regel) schnell und bequem zu 
berechnen, indem wir sie in andere Integrale mit endlich langem Integrations- 


1) Institut fiir Praktische Mathematik der Technischen Hochschule. 
*) I. SzaBé, Die achsensymmetrisch belastete dicke Kreisplatte auf elastischer Unterlage, Ing.- 


Arch, 19, 128 (1951); Beitrdge zur Theorie der achsensymmetrisch belasteten dicken Kreisplatte ins- 
besondere bei elastischer Lagerung, Ing.-Arch. 19, 348 (1951), 
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_ weg und regularem Integranden tiberfiihren. Wir wollen dabei aber die Aufgabe 
verallgemeinern, indem wir folgende Integrale betrachten: 


a 2 

Onl a) = f 2a at mit a=+|al, (1) 
; 

R,,(+) -/{~eia mit a=+|al, (2a) 
0 

ae = \ a | aS gomit J.(t,—<0rund 2,20. © (2b) 


Der Index m soll positiv ganz sein. Die Integrale Q,,(—a) sind als Cauchysche 
Hauptwerte zu verstehen. 


2. Herleitung einer neuen Integraldarstellung fiir Q,,(+ a) 


Eine funktionentheoretische Auswertung der Integrale (1) wird dadurch 
erschwert, dass im Integranden das Quadrat einer Bessel-Funktion auftritt. 
Fiir eine unmittelbare numerische Integration ist ausserdem der unendlich 
lange Integrationsweg lastig. Um beide Ubelstande zu beseitigen, ziehen wir 
die bekannte Integraldarstellung 

2/2 
P(t) = (-1)" 2 fe cos (2 m B) Jy(2 t cos) dd (3) 


0 


heran und erhalten damit fiir Q(-_ a) Doppelintegrale. Vertauschung der 
Integrationsreihenfolge ist ohne weiteres méglich, wenn vor a das positive Vor- 
zeichen steht, da dann der Integrand im gesamten Integrationsgebiet regular 
ist und fiir ¢ > co geniigend stark gegen Null strebt. Es ist also 


7/2 


Onl+ a) = (-1)" = | cos (2 md) [eS ae. (4) 


0 0 
Die Integration iiber ¢ lasst sich jetzt ausfiihren und fiihrt zu dem Ergebnis 


m/2 
Om(-+ a) = (-1)" i cos (2 m 8) [So(2 a. cos) — Ny(2 a cosd)] dd. 


0 


Darin ist S, die Struvesche und N, die Neumannsche Funktion nullter Ordnung. 
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Aus der bekannten Beziehung 


n/2 


J,A2 a COS 7) cosy o dd == =. i 24) Jie t y2(4) 


0 
folgt auf Grund von 


1 0 7 
Ng = on N on Wu "— (=el)" Jiu ae 


die hier fiir uns wichtige Beziehung!) 


[ Nol2.a cos) cos(2m 9) dh =F [Jnu(a) Ng (a) + F-m(4) Ny (a)] | 

0 (5) 
= (-1)" % Ju(a) Na) | 

Damit erhalt man schliesslich 

[. Ini e <= Tela) Nala ym | is. 2a cos) cos(2m 3) do. (6) 


0 


Damit haben wir fiir Q,,(+ a) eine Darstellung gefunden, die fiir die numerische 
Auswertung wesentlich angenehmer ist als (1). 

Bei dem Versuch, fiir Q,,(—a) eine ahnliche Darstellung zu finden, stért 
zunichst der jetzt im Integranden auftretende Pol ¢ = a. Da der Integrand, ab- 
gesehen von diesem Pol, als eine in der komplexen ¢-Ebene regular analytische 
Funktion aufgefasst werden darf, lasst sich der Cauchysche Hauptwert (1) wie 
folgt darstellen 


On(-a) == f 22 piso odie + [feta 


0 
a oe send 
Das Zeichen ~V~ bzw. —\_ soll bedeuten, dass der Weg unterhalb bzw. ober- 
halb am Pol t= a vorbei von Null nach Unendlich gefiihrt werden soll. Auf 
diesen beiden Wegen ist der Integrand regular analytisch, so dass wir jetzt 
dieselben Uberlegungen wie bei Q(-+ a) anwenden diirfen. Man erhalt dann 
mit der Integraldarstellung (3) 
7/2 


Oued) = (ye tf easy | fetes (2 tcos 8) 5, os Jul2 #0089) 4] a) 


0 0 
eS se 
(7a) 


1) J. Dorr, Zwei Integralgleichungen erster A rt, die sich mit Hilfe Mathieuscher Funktionen 
losen lassen, ZAMP 6, 427 (1952). 
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Nun gilt 
1 1 1 9 
Jolt) = = HRY) + > HERD). 


Aus dem asymptotischen Verhalten der Hankel-Funktionen folgt dann 


“Wh Ze 0) : (1) ( 
ee ae di = 24 H\)(2 acosd) 4 | eae 


t—a a 


Cox 


Dabei denken wir uns H}))(k t) bzw. H'?(k t) eindeutig gemacht durch einen von 
t = 0 nach ¢ = —oo iiber die reelle Achse gelegten Schnitt. Beim ersten Integral 
auf der rechten Seite der letzten Gleichung liegt der Integrationsweg auf dem 
oberen und beim zweiten Integral auf dem unteren Schnittufer. Setzen wir 


auf der rechten Seite in der letzten Gleichung ¢ = —z und beachten noch die 
Umlaufrelationen fiir Hankel-Funktionen, so folgt schliesslich 
(2 ¢ cos #) barf i p H\(2 2 cos6) 
re = 71 H\)(2 acos8) ih ee - dz 
0 
1 f HY...) 
inka eager e ay) 
0 
: ao 22c0s# 
= 21 H(2 acosd) — | ae 
0 
Mit denselben Uberlegungen erhalten wir 
Ae CSA —2 1 H®)(2 a cosd) pS" az. (7c) 
t—a z+ a 
0 
Mit diesen semen nimmt (7a) folgende Form an 
7/2 foe) 
(2 z cos B) 
Q,,(—a) = —(—1) fee (2 m #) be (2 acos®?) + 3 Aeros | dd . 
0 


Daraus folgt aber mit den bereits hergeleiteten Beziehungen (4), (5) und (6) fiir 
beliebige Werte a = + || 


0/2 


f Ta gg # 7.(a) Np(a) — (—1)” i? S,(2.@.cos#) cos(2m 3) d®. (8) 


pang im Z 
0 
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Durch Addition und Subtraktion von (6) und (8) erhalt man noch die Bezie- 


hungen: 


Fantini. 
0 
Sa 2 —_4)m ei: 
ee jet ae | So(2 a cos#) cos(2 m #) dd. (9b) 
0 


Auch diese beiden Beziehungen gelten fiir beliebige Werte a = + |a|. Die 
beiden links stehenden Integrale sind wieder als Cauchysche Hauptwerte zu 
verstehen. Fiir m = 0 muss a + 0 sein. 


3. Herleitung neuer Integraldarstellungen fiir R,,(+a) und R»p(-tn,m) 


Der Integrand von Q,,(+ a) ist auf dem gesamten Integrationsweg regular 
analytisch. Deshalb gilt 


af JR 
0 


Daraus folgt unmittelbar 


i! ToT dt = aa = [Jn (2) Nm—1(@) ngs Jin +1(@) N,n(4)] | 
4 7/2 , (10) 
— (-1)™ 2 f E — S,(2a cos) | cos (2 m #) cost dd. | 


Zur Berechnung von R,,(—ty,,) betrachten wir zunachst den Ausdruck 


fo, 3) 


1 ofS) AC, cin fre ; 
i Ga) dt = pet Or ser dé mit a=+la|. 
0 


ams = 

Bei diesen Integralen wird also der Integrationsweg unterhalb am Pol ¢ = a 
vorbei von t = 0 nach ¢ = oo gefiihrt, so dass der Integrand auf dem gesamten 
Integrationsweg regular analytisch ist. Mit Hilfe der Integraldarstellung (3) 
und Beziehung (7b) erhalt man dann 


fe au ieee a . 
| i or sidheaaa Nae! 8 lo =| cos(2m #) |i Hi(2 acos8) 
0 0 
Oe 
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eo : H) 
Teese {ar U Jn (@) m (a) *- Onn a)} 


= + Rn(+@) +04 {Jmn(@) Hi a(@) — Jmsr(a) HY)(a)} 


Unter der Voraussetzung a =1?,, », gilt 


[2 
oe Cin2 |S — $4(2 tp, m cos 9)| cos (2 m 8) cos dd 
0 


5 a Tmsalen, m) Nate, BY : 


Damit haben wir auch fiir R,,(—t,, ,) und R,,(+ a) Integraldarstellungen ge- 
wonnen, die eine mithelose numerische Auswertung erlauben. 


Schlussbemerkung 


Die Struveschen Funktionen S,(x) und S,(x) sind zur Zeit siebenstellig ver- 
tafelt!) fiif den Bereich x = 0(0,02)16. Dieser Wertebereich ist fiir die in der 
Elastizitatstheorie auftretenden Problemstellungen nicht umfassend genug. 
Die Berechnung von S,(x) und S,(x) = 2/a — So(x) fiir gréssere x-Werte ist 
aber relativ bequem mit Hilfe der asymptotischen Entwicklung 


2 12 12-32 12-32-52 
Sola) ~ Nol) + [t-te et | 
_ moglich. 
Summary 


In the mathematical analysis of elastically supported, thick circular plates, 
certain improper integrals occur which must be evaluated between the limits 0 
and oo. The integrals consist of the square of a Bessel function multiplied by a 
fractional rational function. These integrals are transformed into other integrals 
having a finit integration path and regular integrands, in which Struve functions 
of the first and second order replace the squared Bessel functions. 


(Eingegangen: 17. September 1952). < 


1) G.N. Watson, A Treatise on the Theory of Bessel Functions, 2. Aufl. (Cambridge University 
Press, London 1944), S. 666. 
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Le calcul de flutter en régime supersonique 


Par WERNER RupPEL, Madrid, et ROBERT WEBER, La Haye 


Introduction 


L’établissement de l’équation du flutter, pour le calcul de la vitesse critique 
de flutter d’une surface portante d’un avion, nécessite la détermination de 
trois genres de coefficients: les coefficients massiques, les coefficients élastiques 
et les coefficients aérodynamiques. La détermination numérique de ces coeffi- 
cients est assez longue et délicate, a part celle des coefficients élastiques qui, en 
général, peuvent se déduire directement des coefficients massiques. Le calcul 
des coefficients massiques ayant déja été traité ailleurs {1]1), la présente étude 
se limitera a la seule détermination numérique des coefficients aérodynamiques. 

Les paramétres aérodynamiques apparaissant dans les formules pour les 
coefficients aérodynamiques sont des fonctions transcendentes de plusieurs 
variables et complexes. Généralement le calcul de ces coefficients se fait a l’aide 
de tables numériques pour les paramétres aérodynamiques [2], et le volume 
de ces calculs est assez grand. Pour pouvoir cependant les effectuer dans un 
temps raisonnable, on a souvent recours a des méthodes mécanographiques. 
Mais, comme l’usage des tables numériques présente certains inconvénients pour 
la détermination mécanographique des coefficients aérodynamiques, il est néces- 
saire de developper une méthode de calcul évitant l’emploi de ces tables. 

Le présent rapport a pour but de donner une telle méthode de calcul mécano- 
graphique des coefficients aérodynamiques. La méthode est exposée pour le cas 
général d’une aile en fléche munie d’un aileron indéformable en lui-méme. 

On notera cependant que la méthode n’est exacte que pour le cas d’une aile 
plane et pour certains cas limites de l’aile tridimensionnelle qui, dans le super- 
sonique, peuvent étre traités rigoureusement au moyen des paramétres pour 
l’aile plane (par exemple l’aile en delta). Pour le cas général d’une aile tridimen- 
sionnelle, les paramétres aérodynamiques ne sont en effet pas connus, et on est 
alors obligé de recourir 4 la méme méthode d’approximation que dans le sub- 
sonique: a l’intégration par tranches planes successives, en négligeant ainsi 
l’effet du bord extérieur libre de la surface portante. Cette approximation semble 
donner, comme dans le subsonique, des résultats assez satisfaisants. 

La méthode qui va étre exposée et qui est indépendante de tout systéme de 
notations, présente l’avantage que la fréquence réduite, apparaissant dans les 
formules pour les paramétres aérodynamiques comme variable auxiliaire, peut 
étre complétement ¢liminée. I] ne reste plus, de ce fait, dans l’équation du flutter 


1) Les chiffres entre crochets renvoient a la Bibliographie, page 144. 
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que les seules inconnues intéressantes pour le probléme: la fréquence de flutter et 
le nombre de Mach critique correspondant a la vitesse critique de flutter cherchée. 
De plus, comme on peut éviter l’emploi de tables numériques, on évitera égale- 
ment les erreurs numériques d’interpolation qui se présentent toujours lorsqu’on 
utilise de telles tables. Enfin, la méthode est applicable a toute surface élastique 
exposée a un écoulement supersonique (par exemple une pale). 

Remarquons encore que la méthode pourra étre étendue au régime sub- 
sonique lorsqu’on connaitra un développement des paramétres aérodynamiques 
dans ce régime ; il sera toutefois nécessaire alors d’isoler la singularité a l’origine. 


1. Dénominations 


Abréviation Dimension | Désignation 
z m | Coordonnée dans le sens de l’envergure 
Higa m Largeur d’une tranche 
U(z m Profondeur totale du profil 
#1 m Coordonnée dans le sens de la profondeur 
Hol m Ordonnée sur l’axe des ¥/ du bord d’attaque du profil 
tb m Profondeur de l’aileron 
el m Distance du bord d’attaque de laileron a son axe de 
rotation 
a) = Angle compris entre l’axe des » et la direction d’une 
nervure 
Oo = Angle compris entre l’axe des z et le bord d’attaque du 
profil 
Om = Angle de fléche moyen du profil (valeur moyenne des 
angles du bord d’attaque et du bord de fuite avec 
l’axe des 2) 
[(#, 2) gems s2 Répartition de la masse du profil par unité de surface 
Y,(%, 2) _m i-ieme déformation propre de la surface portante 
Yo(2) m Amplitude de translation du bord d’attaque 
B(z) = Angle de rotation des nervures par rapport au plan ¥ ¢ 
YR(2) = Angle de rotation de |’aileron par rapport au plan des 
nervures 
IT, (%, 2) kg m~* Pression aérodynamique exercée en un point de la 
surface portante par suite de la déformation Y;,(¥%, 2) 
de la surface 
y s_* Fréquence de flutter 
v; Sze Fréquence propre de la i-iéme déformation Y;,(%, 2) 
00 lke yoni BC Masse spécifique de lair 
V ms-t Vitesse du vent (dans la direction de l’axe des *) 
c mst Vitesse du son . 
= Nombre de Mach 
wo, =v1/(2V) = Fréquence réduite (définition KUssNER) 
t s Coordonnée du temps 


LO ee ee eee ee 
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2. Equation du flutter 
2.1. Déformation 


Considérons une section d’aile paralléle 4 la direction des nervures (voir 
figure 1). Nous supposerons les nervures indéformables en elles-mémes lorsque 
Vaile effectue un mouvement vibratoire. L’amplitude Y d’un point P(z, x) 
quelconque de cette section sera alors donnée par |’équation [3]: 


Y(z, x) = Yo (2 + A) + B(z+h). (1) 


5 


En développant yp (z + A) et B(z + A) en une série de TAYLOR autour du point 
h = 0, et en remplagant / par sa valeur: 


Cos by Sin 6 


h=l(t— x) <a — ah 


qui s’obtient a partir du triangle PQR de la figure 1, l’équation (1) pour 
Y(z, x) devient: 


Veda z Sree Pel =F Ln Ze) - xk. @ 


Cette équation représente la déformation d’une section d’aile paralléle a l’axe 
des x. 

Les coefficients y,(z) sont sans dimension et s’expriment a l’aide des données 
initiales. On a: 


alt) = [B(2) + 96(2) sind@)] Se, 


lpr " sin 6(z) ] cos? d9(z) sin d(z) 
Y2(2) = | (z) + Vo (2) 2 | c0s*[3, (2) = 6(z)] , 


1 pn m ind z 
yale) = [> A") + 62) = = oar aon pee 
avec 
¥o(2) s- ib ’ ’ etc ’ 


et, d’une facon générale, 


1 d’—1B(z) d"yo(z) | cos’ J9(z) sin’-1 6(z) (r=1) 


Yr(2) = le = 1)l, Saaeeee Re dz cos"[d9(z) — d(z)] 
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Le cas particulier de l’aile droite s’obtient en posant 6, = 0 et celui des 
nervures paralléles a l’axe des x en posant 6 = 0. 

Dans ce qui suit, nous nous limiterons 4 un nombre fini de termes de la 
série (2) et considérerons de plus une aile comprenant un aileron effectuant un 
mouvement de rotation d’amplitude yz autour de sa charniére. Nous suppose- 


“a 


if 


Fig. 1 Fig. 2 


rons l’aileron indéformable. La déformation d’une section d’aile paralléle a 
Vaxe des x s’écrira alors (voir figure 2): 


Ro1 
Yee.s) =| Sm) 2 (x ae +2 le HL bt 8) (3) 
r=0 
{ 0 pour OS fa = 78 
avec 4 = | 


| pour Dr Sail 


Dans le cas ot l’aileron reste toujours horizontal, on posera: 


Yr=—M—4(L—t+e) 2. 


2.2. Les coefficients de l équation du flutter 


On sait que la vitesse critique de flutter s’obtient par la résolution d’une 
équation donnée sous la forme d’un déterminant complexe: 


Pin — Mix + Ain = 0. (4) 


s ° s L : LR oy , ieee ae 
Les coefficients de cette équation sont généralement définis par les formules 


suivantes: 
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Coefficient massique 


| [ [te Y;,(x, 2) w(x, 2) Y,(x,2)dxdz pourt=k, 


0 pour? + k; 
Coefficient élastique 
V2 Mii pour? =k, 
Ve — : 
| 0 pour? + k; 


Coefficient aérodynamique 
Ape I / (2) Y,(x, 2) I(x, 2) dx dz, 


Y,(x, z) et Y,(x, z) représentant les formes propres calculées de la surface por- 
tante ou encore les formes de résonance obtenues dans un essai au sol. Les inté- 
grales dans ces formules s’étendent a toute la surface portante. 

On voit que pour calculer le coefficient A,, il est nécessaire de connaitre la 
pression J7;, exercée en un point de l’aile par suite de la déformation Y,, de 
l’aile. Or cette pression JJ, n’est connue explicitement que pour le cas d’une 
aile «plane», c’est-a-dire une aile a bords paralléles et d’envergure infinie. Pour 
un cas pratique de calcul de la vitesse critique de flutter d’une aile d’envergure 
finie, on est alors obligé de recourir 4 une méthode d’approximation qui con- 
siste a diviser l’aile en un certain nombre de tranches par des sections paralléles 
ala direction du vent. On suppose que ces tranches d’ailes ont leurs bords paral- 
léles et que la déformation est constante dans chaque tranche. Cette méthode 
conduit a la définition suivante du coefficient A; ,: 


1 


Aj = = ys In(en) | Y,(%, Zn) IT, (x, Zn) i Az, , (5) 
0 


n 


1,(z,) désignant la profondeur moyenne de la tranche n, Y;(x, z,) désignant la 
déformation «moyenne» de cette tranche et J//,(x, z,) désignant la pression 
exercée en un point x de cette tranche par suite de la déformation Y;,(x, z,). 

Dans ce qui suit, nous nous occuperons uniquement du calcul du coefficient 
A», le calcul des u;;, et ~;, pouvant se faire de la maniére habituelle, c’est-a- 
dire en divisant l’aile en tranches paralléles aux nervures. 


, 


2.3. Calcul du coefficient aérodynamique A, , 


D’aprés l’€quation (3), on voit que la déformation Y, peut s’écrire sous la 
forme d’une combinaison linéaire des termes y;, qui s’obtiennent a partir des 
données initiales du probléme. On montre que la pression I7;,(x, z,) peut égale- 
ment s’écrire comme combinaison linéaire des termes y;,, correspondant a la 
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déformation Y,, (voir Annexe, page 137): 
R 
IT, (x, Zn) = Valen) p;(*) (6) 


En introduisant les équations (3) et (6) dans l’équation (5), on obtient pour A;,: 


ta Bef E32 wa) (Sn) 


agi =o Vin ft = Xo = Law — 6) 1S, ra bt 


1-t 


et, en transformant légérement cette équation: 


R= oR 
Ai, = ye EFS 7 Yes n fig ae Xo)" ps dx 
. (7) 
aes 
+ 2S vinverle | 2e- Rp La tie) pax Az, : 
1-t 
Les intégrales apparaissant dans l’équation (7) peuvent étre ramenées au 


cas particulier de l’aile droite (6,, = 0), c’est-a-dire d’une aile ayant ses bords 
perpendiculaires a la direction du vent, par la relation: 


Ps is COS On, (Ps) 6m =0 2 


En introduisant encore les abréviations suivantes: 


af 
4 1 
ly | 2 (# — 0)" (Ps)on=0 4% = 00 V? Tira 2 ORT SR-1) a) 
0 
1 ' ; 
by f 2 (e— xo A+ 42) Drlay-0 dt = CV? eZ tne, (BD) 


{1 
ou les paramétres m,, sont des paramétres aérodynamiques qui peuvent étre 
explicitement calculés (voir Annexe, page 141), l’équation (7) pour A, ;, devient: 


Aer = (L)* 0 :088n ot S (2) (3) 3 emer} ar 


ju? ¥=O'S 


ou bien, en notation matricielle: 


Ain = (4) 00 085m pez & (2) Weel} 4e- 0) 


V2 ie 
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Dans cette équation: 

(m,,] est une matrice carrée, représentant les paramétres aérodynamiques de 
Vaile pour la tranche ; 

[y;,] est un vecteur contrevariant (ligne) représentant la 7-iéme déformation 
de l’aile pour la tranche n; 

[ysx] est un vecteur covariant (colonne) représentant la k-iéme déformation 
de l’aile pour la tranche n. 

Lorsque l’aile 4 étudier a des nervures paralléles a l’axe des x (6 = 0), 
l’ordre de chacune des trois matrices ci-dessus diminue d’une unite. 

Les paramétres aérodynamiques m,, sont représentés par des fonctions con- 
tinues de la fréquence réduite w, (du moins pour des valeurs finies de w,). On 
peut donc développer ces fonctions en des séries convergentes autour du point 
w, = 0. Ces développements sont connus sous la forme générale suivante (voir 
Annexe, page 143): 


-5 mM, s, o cee | (,) (10) 


7 étant l’unité imaginaire. 

Les coefficients m,,,, sont des fonctions du nombre de Mach et, pour les 
parameétres de l’aileron, de certaines grandeurs géométriques de l’aileron. 

Le nombre de termes nécessaires de cette série dépend d’une part de la 
précision désirée et d’autre part de M et de w,. Pour les problémes stationnaires 
(divergence, inversion, stabilité, etc. [4]) caractérisés par des valeurs de w, extré- 
mement faibles, ainsi que pour les problémes de flutter a trés haute vitesse, 
on pourra se limiter aux deux premiers termes de la série. Pour les calculs de 
flutter usuels (1 < M S 2) il faudra en général prendre les huit premiers 
termes de la série (0 So <7). 

En introduisant la série 4) dans l’équation (9) on obtient: 


Ay, = (;)" 00 COSOn, Tat 1 d dy | —J @,) "(2 =)" [vir] ™rs,0] ¥en] Aen 


ou bien: 
Ay, = (~)’ 0o COS Oy, 7 1 | 
yf 2 On \oo by NO ay) 
«FE (A (EY bellmacllread dee | 


In inversant la suite des sommations dans cette équation et en remarquant que 


le rapport 2 w,/l, = _v/V ne dépend pas de la tranche n considérée, on a finale- 
ment pour A,,: 


Ash = 09 C088 eG. i) “(+ aa A, (12) 


a 
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avec by\o +2 
oo (2) Bo, Az (13a) 
et 
Bon = [Vir] [Mrs, ol lyse] (13b) 


ou, en séparant en parties réelle et imaginaire: 


/ ® i 
Ay, = Ain +t ik 


avec 5 4 o [|v \20-2 
Ais = 000088 a= 5) (-1) (7) Ase 
" 4 a yp \2o-1 
Aj, = 20 cos Om /M?—1 pa, ( 1) a (+) Aros 2 


Les indices 1 et o de l’expression (13b) pour B,,, indiquent qu'il faut cal- 
culer B,,, pour chaque o et chaque tranche. D’aprés l’équation (13a), la somme 
de toutes les tranches donne l’expression A, qui ne dépend donc plus que de o. 
L’équation (12) donne ensuite le coefficient A;,, en fonction du rapport v/lV. 


Remarquons que 
matricielle suivante: 


Ajx, = Q9 60S 


’équation (10) peut également s’écrire dans la forme 


v 


V 


In 


ails 


(14) 


By Bog] [429] 


beater (i 


[(- Us (F)] étant une matrice 4 une ligne; 

(ey BF étant une matrice rectangulaire ; 
[Az,,] étant une matrice a une colonne, 
et. B,,, étant défini par l’équation (13b). 


Enfin, en remplacant dans les équations (12) et (14) V par M c, on obtient: 


Lp = Og 6* COSOm, yet Zh (—j) 07-2 A, (15) 

avec = 1, \ar2 = 
Ay SS a (= Be Ae. , (16a) 
53 = [yir] Raa [vse] (16b) 
ou bien ; ee ee | B 
Asx = Qo 64 COSSm Vices Yar | (5%) B,,| [4z,]. (17) 


On voit qu’en employant pour les coefficients A;, les expressions (15) ou 
(17), ’équation du flutter ne contient plus que les deux inconnues » et M. 
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3. Résolution de l’équation du flutter 


Les coefficients 4; 4, Pix et A; une fois calculés, il ne reste plus qu’a résoudre 
l’équation du flutter (4). Deux méthodes sont possibles, selon que l’on choisit 
comme équation de définition de A,,, une des équations (12) ou (14) contenant 
y, V et M comme inconnues, ou bien une des équations (15) ou (17) ne conte- 
nant que v et M comme inconnues. Par suite du caractére compliqué des coef- 


Meo/ 


subsonigue SUPETSOIN GUE 


0 85 1 f7 2M 3 4 M,, 


ficients A,;,, on sera, dans les deux cas, obligé de recourir 4 un procédé de réso- 
lution semi-graphique de l’équation (4). 

Dans le premier cas [équation de définition (12) ou (14) pour A;,], on con-- 
sidérera M comme paramétre et » et »/V comme inconnues. On se donnera 
différentes valeurs de M et on déterminera, pour chacune de ces valeurs, les 
inconnues v et V de telle fagon que l’équation (4) soit satisfaite, la solution du 
probléme étant celle pour laquelle la valeur calculée de V correspond a celle de 
M introduite comme paramétre. On part donc d’une vitesse critique supposée 


Veuy et on détermine d’abord le nombre de Mach effectif M correspondant, par 
la relation: 


COS 0» 
M = Viun rere 
Veup devant étre choisi tel que M soit supérieur a 1. 
Les deux équations résultant de l’équation complexe (4) ne contiennent alors 
plus que » et »/V comme inconnues. Pour la résolution de ce systéme d’équa- 
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tions, on procédera comme dans l’incompressible, c’est-a-dire on calculera les 
coefficients des équations pour un certain nombre de valeurs »/V, et on déter- 
minera y graphiquement de telle sorte que les deux équations soient satisfaites. 
On obtiendra ainsi les solutions y et »/V de l’équation (4) correspondant a la 
valeur de M introduite dans le calcul. Du couple de valeurs calculées y et »/V, 
on déduira V que l’on appellera vitesse critique calculée V,,;. On effectuera 
ainsi plusieurs calculs, chaque fois en partant d’une autre valeur pour V,,», de 
fagon 4 pouvoir tracer la courbe V.,, en fonction de V,,,,. L’intersection de la 
droite Vig, = Vy, avec cette courbe donnera la vitesse critique de flutter 
cherchée (voir figure 3). Si la courbe posséde plusieurs points d’intersection 
avec la droite, l’intersection correspondant 4 la vitesse la plus faible donnera 
la vitesse critique cherchée. 

Dans le deuxiéme cas [équation de définition (15) ou (17) pour A;,], on cher- 
chera pour différentes valeurs de M la ou les valeurs réelles », annulant la 
partie réelle de l’€quation (4) et la ou les valeurs réelles », annulant la partie 
imaginaire de l’équation (4). On tracera ensuite les courbes », et v, en fonction 
de M. L’intersection des deux courbes donnera la fréquence critique ¥,,,;4 et 
le nombre de Mach critique M,,,, correspondant a la vitesse critique de flutter 
cherchée; la vitesse critique elle-méme s’obtiendra par la formule 

V. c 


crit _ *orit ve 
COS), 


ANNEXE 
Détermination des paramétres aérodynamiques 


On sait qu’en régime supersonique les effets ne se répercutent que dans le sens 
de l’écoulement. II s’en suit que la pression aérodynamique JI(*) exercée en un 
point x d’une aile plane A, ayant une profondeur égale a t/ et effectuant un 
mouvement d’amplitude Y, est égale a la pression exercée en un point (¥ + 1— 7) 
de l’aileron d’une aile plane B dont la partie située en avant de l’aileron n’effectue 
pas de mouvement, l’aileron en question ayant une profondeur égale a t/ et 
effectuant le méme mouvement d’amplitude Y que l’aile 4. Il suffira donc de 
calculer les paramétres m,, pour une aile plane de profondeur Tt 1, tous les para- 
métres m,, pour une aile plane de profondeur totale 1 et munie d un aileron de 
profondeur t/ pouvant étre déduits directement des parameétres ™,.. 

De plus, comme les paramétres aérodynamiques pour une aile plane dont les 
bords font l’angle 6,, avec la direction du vent se déduisent de ceux d’une aile 
droite (6,, = 0) par simple multiplication par cos dn, On peut se limiter, pour la 
détermination des paramétres ,,, au Cas particulier 6,, = 0; on pourra par consé- 

lement supposer %» = 0. 
ar eae onc une aile plane droite ayant une profondeur totale 
égale a t/ et effectuant un mouvement vibratoire, harmonique en fonction du 


temps, d’amplitude: 


_— 
art 
~~ 


R 
par ih ; 
ites bee rT yt pivt 
Y(%) e7%* = : », Di awe ees 
T= 


138 WERNER RupPEL et ROBERT WEBER ZAMP 


cette aile étant exposée a un écoulement paralléle animé d’une vitesse super- 
sonique V. 

Pour les besoins du calcul de flutter d’une aile il sera suffisant de pousser la 
série (1) jusqu’au troisi¢me terme (R = 2). 

D’aprés la définition des paramétres m,, et d’aprés ce qui a été dit ci-dessus, 
il nous faudra calculer des intégrales de la forme générale suivante: 


NA CHE if 2"(~ — &)" I(x) dx , (Oly SZ) a) 
0 
€ étant un point de référence quelconque. 


1. Calcul de Mo(t) 


La pression J/(t) exercée au point t de la tranche d’aile est donnée par la 
formule [5]: 


(zt) = 2 0, V (Meni | 
2 M2 ; 1 2M? (3a) 
x - w(t) + Mi 1 o, | w(x) i BT wanes O, (t = x) ts | 
0 
F 2 on = = 
sh Serr oe, [yo Bolt) + v1 Palt) + Ys Palt)] (3b) 
avec ¥ 
t K(v, u) = ve-I*[f,(u u) +7 v Jo(v u)) 
(=) 
i 1 d es dY HR aes 
wht) epee aE lath pier ap raed te 


Jo et J, sont les fonctions de Brssev d’ordre 0 et 1. En posant la série (1) dans 
w(*), et sachant que 


vy 20, 
Caer Bs 
on obtient: 
w(*) = Yo7 Wy, + yy (L +27 WX) + Yq 4 (% + 1 W, ¥*) . (4) 
On sait que l’on a pour M(t): 
M(t) =1 | IT(x) dx (5a) 
0 
aw V2 ee 
a VM?2—1 


2M? is ye 2M? 
xX 1)— w,(t) + Fe a | W,(*) K( 5 Moi ,(T *)) ax 


t | (5b) 
| 


0 


A 4 ine aes i 
= 0) V? Vue = y Zz lve Moo(T) + Vy Mo1(T) 4 Vg Mo(T) | (5c) 
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avec 
x 


2 


w, (x) = ie we) d= rojo tnletio)+n4(F rio, +). (6) 
2 


Le paramétre M(t) est connu et admet un développement en une série conver- 
gente autour de w, = 0 ([2]). On a: 


co fore) 
Migg(t) = 327 ho,o(M) (—] wy) = YS” Moo,c (—7 wp)” (7) 
=0 = 
avec 5 sie 
M2 
Fog — 0, hoya ~ 6 (M?— 1)8 (age 1), 
M4 
aa aeer Won "19 Ta aye Me tS) 
M4 : 
oh eet Vg Seana Ta hoe 60 (MP = is (8 M4 +12 M?2+1), 
M2 Me 
L = — wee 4 2 
fils Cie aye’ Re 180 (ME (8 M4 +12 M?+5), 
Ms 
hos (64 Ms + 240 M4 + 60 M? + 5) 


8 5040 (M2 — 1)? 
ou bien, d’une facon générale: 


Mo-2 (a/2) —1 (2 M)e-2-2m 


: = 
fool M) (¢,=1) (M@? = 1)5-1 x Ga =o mite ee 


Remarquons ici que cette série est également valable pour des valeurs complexes 


de o,. 
En comparant les équations (3a), (3b) et (4) avec les équations (5b), (5c) et 


(6), on voit que l’on a: 
Po(t) + 2 moo() (8) 


Pilt) = j @, 


P2(t) = 4 Mox(7) - (9) 
En intégrant l’équation (8) de 0 a t on obtient: 


et 


M(t) = a Moo(t) + 2 [rat 


et, en introduisant pour 799(t) la série (7): 


[o-e) oo) ott ; 
moy(t) = — 372? hoo(M) (—7 @y) a = fio, o(M) (—f ,)%. (10) 
o=0 


_ Sachant que fy o(M) = 9, l’équation (10) peut s’écrire: 
Mg, (7) a a (—7 @,)° eae? ho,o(M) — ho,o41(M)} 
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Nous posons encore: 


2 
hy o(M) mre ey ho, o(M) — ho,o41(M) (11) 
et mo,(t) sera égal a: 
co foe] 
Mg(t) = StF hy, o(M) (—j @,)? = D Mion, o (—7 or)? - (12) 
o=0 o=0 


En intégrant l’équation (9) de 0 aT, on obtient: 


Mo9(T) ae 4 | Wiest) de 
0 


et, en introduisant pour 77%9,(¢) la série (12): 
nt qzot2 
. Mo9(T) = > 4 ae a Ay o(M) (—] @, = Sino | a ,)? (13) 
2. Calcul de M,(t, §) 


s) <1 [2 (x) Ma) ax 
a (0 


et, en effectuant une intégration partielle, 


On a pour M(t, &): 


0 
= op Vit. — 4 __ fp, miele, €) + v4 Hult, 8) + 7, Malte Oe 
2 yM2 on 11 ’ 2 12 ’ 
Il s’en suit: 
x 
My ,(t, &) 


I 
dS 
a 
= 
o 
a 
bho 
et 
S| 
QQ 
ad 


On aura donc, en introduisant dans cette équation tour a tour les séries (7), (12) 
et (13): 


co le. °) 
— 3 Tt . 
Myo(t, 8) = 32 (4 — § — = 2 ho, o(M) (—7 a4)" =D Mio,o (—7 )”, 
o=0 o=0 
foe) 7 lo.) 
malt, €) = 3/2(6—§ — Ta) e048 hy o(M) (—j 0)” = DY Fa (—7 04)", 
o= o=0 


T tote iad 


M(t, €) ay (= g ‘oe ) 4 o+2 hy o(M) (=] (,)° = 2, Mis,0 (—7 Wy)? « 


Vol. IV, 1953 Le calcul de flutter en régime supersonique 141 


3. Calcul de M2(t, &) 
On a pour M,(z, €): 


MGs, 8) = [a (4 — &)? I(x) dx 
0 
et, en effectuant deux intégrations partielles successives: 


M,(z, &) =4 c — 1)? M(t) + 2 (€ —72) | Male) dx +2 [ [ee dt a 
0 0 0 


1 4 ee a ao 
Co ie : yme—1 [Yo Mao(t, €) + Yi Mai(T, €) + Yo MoQ(T, Cals 
Il s’en suit: 
i 4(t, §) = 4 le 9 +2 (8-7) | Higg(x) de +2 | f rigg(t) a ea 
0 0 0 


On aura donc, en introduisant dans cette équation successivement les séries (7), 
fr 2et (13): 


= Dae ; 
Mind, £) = 374 |(e—8)*— (© — 8) + ea | © aol) (7 0) 


o=0 
= Bag 2G * ok eS 
Mins(t, 8) = D4 te = 88-4 egy | Ot el im) 
=) iria,o (=7o,)?, 
o=0 
= LEG 2 ue 
pena?) pind | EL pals MENG Ra oR) aan 
oi o+2 i 
Ih a(M) (-j oy)? 


co 
= >) Wiss,0 (4 Cy)? ‘ 


4. Calcul des paramétres m,, 


Nous pouvons maintenant calculer sans peine les parameétres m,, et Mr, 
définis par les équations (7a) et (7b) du paragraphe 2 du présent rapport. 

Les paramétres m,, concernant l’aile (0 S1,s 2) se déduisent des fe 
métres m,, en posant t = 1 et & = 0 [puisque le point de référence est, d apres 
l’équation de définition (8a) du paragraphe 2, le bord d’attaque]. On a donc: 

M,, = ™M,,(1, 0) . (O°S r,s S 2) 


Ts 
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Les paramétres mp, concernant l’aile (0 < s < 2) et ayant comme point de 
référence l’axe de rotation de l’aileron s’obtiennent a partir des paramétres m,., 
en décomposant l’intégrale de définition (8b) du paragraphe 2 de la fagon sui- 
vante (avec ¥%) = 0): 


1 1 1 
j [2(x-14+1-2) pydx=l 2 tx—(— a) bax —20f peda 
‘err 1-7 1l-rt 


On a donc pour mp,: 
mp, = M;,(1, 1 —t) —™,,(1 —7, 1 — 1) — 2 [M%o,(1) — M,(1—7)]. (0s 2) 


Pour calculer les paramétres m,p concernant l’aileron (0 Sv S 2 et r= R), — 
nous devons effectuer une transformation des paramétres du fait que la rotation 
de l’aileron ne s’efféctue pas autour de son bord d’attaque, mais autour de sa 
charniére. En désignant par %7,p le paramétre de l’aileron, pris par rapport au 
point de référence &, lorsque l’aileron effectue une rotation autour de sa charniére, 
on a d’une facon générale: 


M,R(t,&) = M,1(t,&) —2em, (7,6). (OS rS2etr=R) 


Pour les trois premiers de ces paramétres (0 < y S 2), le point de référence est, 
d’aprés leur équation de définition (8a) (paragraphe 2), le bord d’attaque de 
Vaile, c’est-a-dire situé & une distance égale a |r — 1| du bord d’attaque de I’aile- 
ron. On a donc & = t — 1, et on obtient: 7 


MR =M,R(t,t — 1) =mM,,(t,t —1) —2em,.(t,t-—1). (0S 7S 2) 


Enfin le quatrieme paramétre mpr, qui correspond au moment par rapport a 
l’axe de rotation de l’aileron, lorsque l’aileron effectue une rotation autour de son 
axe [d’aprés l’équation de définition (8b) du paragraphe 2}, s’obtient de m, p(t, §) 
en y posant = «; on a: 


MRR = M,R(T, €) = My,(T, &) — 2 € M,(T, €) . 


5. Cas d’un aileron avec fente fermée 


Dans ce qui précéde nous n’avons considéré qu’un aileron avec fente ouverte 
(définition KUSSNER-SCHWARTZ [6}), c’est-a-dire avec une compensation aéro- 
dynamique. Il nous faut done encore considérer le cas d’un aileron avec fente 
fermée. On sait que cette caractéristique n’a d’influence que sur le paramétre 
aérodynamique provenant d’une translation de l’aileron, donc le paramétre Woo. 
Soit Mo) le paramétre de l’aileron avec fente fermée lorsque celui-ci effectue une 
translation. Considérons la figure 4. Le profil d’aile A, est rigide et effectue un 
mouvement vibratoire de rotation autour du bord d’attaque ¥= 0. Le profil A, 
se compose d’une suite d’ailerons 4 fente fermée et disposés en escalier; chaque 
aileron a la méme profondeur Ax; les ailerons effectuent chacun un mouvement 
vibratoire de translation, tous avec la méme. fréquence et en phase, de sorte que 
lorsque les ailerons passent par zero (axe des x) le profil A, est une droite. En 
désignant par f l’angle de rotation du profil 4,, on voit que l’amplitude de vibra- 
tion d’un aileron par rapport a l’aileron précédent est égale a 1 B Ax. La force 
aérodynamique agissant sur le profil A, est égale (4 un facteur prés) a 7 B my,/2, 
tandis que celle agissant sur le profil A, est égale a 1] Y B Ax 7Hgo. A la limite 
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WA Ax +> 0, les deux forces aéro- 
dynamiques sont égales: 


i 
Mo, = 2 | a ax. 
0 


Par le changement de variable 
#% = 1 —T, on obtient: 


ul 


Mo, = 2 | a ORS 


0 


Nous développons My) en une 


Fig. 4 série semblable a la série (7): 
Moo a is ho, o(M) (=7 Wy)? o 
o=0 


En introduisant pour m,, la série (12), avec t = 1, on obtient: 


as ae 
Inyo M) = 2 To,o(M) | 29 de = To, o(M) | 


On voit donc que pour un aileron non compensé aérodynamiquement (fente 
fermée) il suffit de remplacer hy ,(M) par 


= o+1 
ho, o(M) = ie, aie hy, o(M) 


dans les formules pour m,p (0 Sr 2etr= k). 


6. Les paramétres m,, pour r,s > 2 


Dans certains cas (pales ou autres surfaces de ce genre), il sera nécessaire de 
considérer plus que trois termes dans le développement en serie de TayLor de Y. 
Il sera alors également nécessaire de connaitre les paramétres m,, pour 1, s > 2. 
Pour déterminer ces paramétres, on posera, en généralisant les formules obtenues 
dans les paragraphes précédents: 

le.e) 
My g = J) mys, (—] @,)° . 
o=0 
Comme on pourra facilement s’en rendre compte en raisonnant par récurrence, 
les termes m,, g S obtiennent par la formule générale suivante: 


2" (o + s) 
My so = aay ee Mg 5,0 (v, s+ R) 
avec 
[f.2 pour s=0, 


hy pour s2l. 


Mos,0= 4 28-18! (0 +1)! 
| (o+ s)! 
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7. Conclusion 


Nous avons ainsi obtenu un développement en série autour de w, = 0 pour les 
seize paramétres m,, entrant en ligne de compte pour un calcul de flutter d’une 
aile avec aileron. Comme on a pti s’en rendre compte, tous les paramétres se dé- 
duisent d’un seul d’entre eux, en l’occurrence le paramétre 7799. De plus, les coeffi- 
cients des séries sont des produits d’une fonction de w, seul et d’une fonction de 
M seul; d’ailleurs cette fonction de M ne doit étre calculée que pour un seul para- 
métre; en effet, il n’apparait dans les coefficients que les fonctions ho, et fy,4 et 
les deux fonctions sont liées par une relation trés simple [équation (11)}. Dans les 
cas pratiques, on déterminera donc uniquement hy g, et on en déduira h,,, par la 
relation (11). 

Tous Jes résultats obtenus dans cette annexe pour les paramétres m,, sont 
rassemblés dans le tableau ci-contre. 
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Zusammenfassung 


In vorliegendem Bericht wird eine Methode entwickelt zur numerischen 
Berechnung der im Uberschallgebiet geltenden Luftkraftkoeffizienten mit Hilfe 
von Lochkartenmaschinen oder programmgesteuerten Maschinen. Dazu werden 
die Ausdriicke fiir diese Koeffizienten in Reihen entwickelt, wobei jedes Reihen- 
glied ein Matrizenprodukt ist. Dadurch wird es ebenfalls méglich, die reduzierte 
Frequenz, die bei Flatterrechnungen immer als Hilfsvariable auftritt, zu elimi- 
nieren. Die Flattergleichung enthalt dann als Unbekannte nur noch die wirklich 
interessierenden Gréssen: die kritische Frequenz (Flatterfrequenz) und die kri- 
tische Machsche Zahl (Flattergeschwindigkeit). 

Da die Luftkraftbeiwerte fiir den allgemeinen Fall einer schwingenden Trag- 
flache noch nicht bekannt sind, wird diese streifenweise eben angenommen. Die 
Methode wird dargelegt an Hand einer schrag zur Windrichtung angestellten 
Tragflache mit ebenfalls schrag zur Windrichtung angebrachten Querrippen. 
Alle anderen Falle lassen sich ohne weiteres aus diesem allgemeinen Fall ableiten. 

In einem Anhang werden die zur Berechnung der Luftkraftkoeffizienten be- 
nétigten Reihenentwicklungen der Luftkraftbeiwerte bestimmt. 


(Regu: le 20 novembre 1952.) 
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Kurze Mitteilungen — Brief Reports - Communications bréves 


Résolution de problémes aux limites au moyen de transformations 
fonctionnelles') 


Par Paut Burcat, Neuchatel 


1. Introduction 


On a utilisé, pour résoudre certains problémes aux limites, la transformation 
finie de FoURIER?). 

Je me suis proposé, dans ce travail, non de prendre pour point de départ une 
transformation et de déterminer son champ d’application quant a la résolution 
des problémes aux limites, mais de rechercher un opérateur ou, ce qui revient au 
méme, un type de série approprié a un probleme donné’). 

La relation de Green, j’en fais usage, joue un rdle fondamental dans plusieurs 
questions relatives aux équations différentielles*). Elle s’écrit (pour les fonctions 
d’une variable) 


(L u,v) — (wu, L* v) = (0, AV = A(U, V), 
si l’on pose 
Lu=p,u™ +p, u®-YV4 +--+ po, u, 


b 


Ciera / fgdx (fet g, fonctions réelles et intégrables) 


a 


et que l’on désigne par L* l’expression différentielle adjointe de L wu et par 
(U, V), le produit intévieur des vecteurs 


U [u(a), u’(a), ..., w"-V(a), u(b), w’(b), ..., wl" D(d)] 
et 


V (v(a), ..., v™-V(a), v(b), ..., v-N()) . 


A est une matrice réguliére a 2 » lignes et 2  colonnes; les fonctions réelles p; 
satisfont aux hypotheses habituelles®). 


2. Probléme 
Considérons l’équation différentielle linéaire du n-iéme ordre 


Ly+ry= f(z) (1) 

1) Résumé des trois premiers chapitres d’une thése élaborée sous la direction de M. Cu. BLanc, 
professeur a l'Université de Lausanne, et présentée, en 1950, A Neuchatel. Voir aussi C. r. Acad. 
Sci. Paris 231, 321-323 et 398-400 (1950). 

2) G. Dorrscu, Math. Ann. 112, 52-68 (1935). H. Knress, Math. Z. 44, 266-292 (1939). Voir 
aussi IDA RoETTINGER, Quart. appl. Math. 4, n° 3, 298-319 (1947). 

3) Cf. R. V. CouRcHILL, Amer. math. Monthly 59, n° 3, 149-155 (1952). 

4) Cf. M. Piconr, Appunti di analisi superiore (Rondinella, Naples, 1940), p. 746. 

5) Cf. E.L. Incr, Ordinary Differential Equations (Dover Publications, New-York, 1944), p.211. 


———s 
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et les conditions aux limites, supposées linéaivement indépendantes, 
L(y) = al? y(a) + aff! y’(a) +--+ + al) y*-Na) 
+ off) y(b) $s + al) yD) = al 


vy et les a sont des constantes réelles; f(%) est une fonction génévalement continue 
(sttickweise stetige Funktion) sur (a, b). 

Nous nous proposons de trouver une série de fonctions propres qui représente, 
dans l’intervalle ouvert (a, b) ou sur le segment (a, 6), la solution y du probléme 
[(1) et (2)], supposé déterminé. Nous n’envisagerons ici que le cas dans lequel le 
probleme homogéne correspondant est auto-adjoint. 


3. Méthode de résolution 


. . eee . = 
Soient v une fonction que l’on va définir et Y le vecteur de composantes 


y(a), ..., y™-V(a), y(b), ..., yb). 
Ona & 
(y, L* 0 + rv) = (f,v) — A(Y, V). (3) 


Dans les questions ov la formule de Green intervient, on cheyche d’ovdinaire a 
déterminer v de maniére que A(Y, V) soit nul; nous choisivons au contraive v de 
facgon que cette expression soit une fonction connue des quantités a), 

M désignant une matrice réguliere quelconque, a 2 ” lignes et 2 ” colonnes, 
posons 


N=(M-2’A: Y=MY, V=NV7; 
nous avons 
A(Y, V) = Y’AV = (M-!Y)/AN-1V=Y"(M-1)AN-V 
= (MA N-“Y, V) = Ey,.V) 1). 


Le second membre de la formule de GREEN prend la forme (dite canonique) 


E’Y,V = 3 %0) 


les Y;,(y) et les V;,(v) étant respectivement les composantes des vecteurs Yet V. 

Jusqu’a présent, la matrice M était quelconque (réguliére) ; désovmais, ses n 
premiéres lignes seront constituées par les coefficients qui figurent dans les conditions 
aux limites (2); nous aurons 


A(Y, V) = EY, V) = Dd ol V, w) + S¥0) ) V;(v) . 


t=1 g=n+1 
Nous prendyons, pour v, une fonction satisfaisant aux conditions 
Vv) = 0 G=n+1,n+ 2,...,20); 


A(Y, V) sera une fonction connue des a) (¢=1,...,”) et la relation (3) deviendra 


(y, Lev +70) = (7,0) — D* a V,(0) . (3° 


1) Y’, par exemple, est la matrice transposée de Y, matrice a une colonne. 
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Pour définir complétement la fonction v, nous l’astreignons a vérifier l’équa- 
tion ; 
L*v+Aav=0 (4) 
outre les conditions trouvées précédemment: 

V;(v) = 0. (7 = 2+ 1.0 52,5. eZ a) 


Nous supposons dés maintenant que le probléme [(4) et (5)], nous l’appellons 
«probléme associé») de [(1) et (2)], est auto-adjoint. Désignons par 4, Ag, ... Ses 
valeurs propres et par 4, V2, ... les solutions propres correspondantes. La relation — 
(3’) nous donne, si nous tenons compte de l’équation (4), 


(Y, x) = or ny al) V;(v,) — (4, Up) |; (3”) 
tI 


(y, U;,) est la transformée de la solution y du probleme [(1) et (2)] par l’opérateur 
intégral de noyau v,. Si les fonctions propres forment, apres avoir été normées, 
un systéme complet?), la formule (3”) donne, a un facteur prés, les coefficients du 
développement de y selon les fonctions orthogonales vj, Ug, .... 


4. Equation du second ordre; conditions de troisiéme espéce 
Nous allons résoudre, par la méthode développée dans le § 3, le probleme: 
Ly ry = (pls) Vl + gla oe ia (6) 
% (a4) + % y’(a)= a, Bs y(b) + Byy’(b)=B, (a} + a8) (B3+ B3) +0; (7) 


p(*) et q(¥) sont des fonctions continues sur (a, 6); on suppose de plus que p(*) 
est positive et qu’elle admet une dérivée premiére continue [sur (a, b)]. 


Gin 
0 pia) 0 o |i 
Pi —p(a) 0 0 0 | 
0 0 0 —p(b) ||" 
| ot Se opin ako, 2m 


Nous pouvons prendre pour M la matrice réguliére 


Oy Xe 0 On 
0 0 Bs By 
|| Oy — 0 0 


0 0 — By Bs 


Calculons la matrice N, puis les composantes du vecteur V; nous obtenons 


Vai-a(0) = (1) atte pla) v(a) + pla) va) , | 
ote 4 (@=1,2) (8) 
V,,(v) = Rt pp (b) v(b) + (—1)# 2 pe PCO) v’(b) . | 


1) (4) et (5) et Lu+Av=0, L,(v) = 0 (i =1,...,) sont des problémes équivalents. Dans 
certains cas cependant (équation du second ordre avec conditions aux limites générales, par exemple), 
le premier conduit aux résultats les plus simples. 

2) Vollstandig (CouRANT et HivBerr); closed. 
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Le probléme associé de [(6) et (7)] est le probleme de StuRM-LiovuviLte. Les 
théorémes d’équiconvergence nous permettent d’affirmer que la solution y de 
[(6) et (7)], fonction 4 variation bornée, est représentée, dans (a, b), par son 
développement a l’aide des v, [on suppose que p(x) et g(x) satisfont aux hypothéses 
de ces théorémes]. Si a, £, + 0, y est égale A son développement sur le segment 
(@, b) (théoréme de Haar). 

Les coefficients a, de ce développement sont données par la formule 


1 1 Ot a , 
ry tr E p(a) (- ane iy Up(@) + eso v(a) 


(9) 
+ B p(b) (a a v;,(b) RB ie 4(0))| = Ay}, | 9 


deduite des relations (3”), pour = 2, et (8); on a posé N, = (v,, 0;,) et 
A, = (1/N,) (f, vz). Les A;, sont les coefficients du développement de la fonction 
{(#), second membre de l’équation (6). 

On peut mettre le second membre de (9) sous une autre forme, plus commode 
dans certains cas. 

Soient g(x, A) et h(x, 4) deux solutions de l’équation Lvu+Av=0 qui 
constituent un systeme fondamental et telles que 


NG a =O AAG, A) =O, 


quel que soit A. Posons 
p(a@, 4) = c(A) et h'(a, 1) = aA) . 
L’équation caractéristique du probléme associé peut s’écrire 


a, (A) oa oy A(A) m1 
Bs g(b, A) + Bag’(0, A) Bg h(b, A) + By h'(b, A) © 


Désignons par @; la valeur commune de ces rapports, pour A= A, . On obtient 


aisément 
i By U,(0) — Bs vx (5) 
On ~~ “(BR + BR [e’(, Ax) A(B, Ay) — B(b, Ay) A’(6, A) 
p() Ba Ue(0) — Bs ve(2) 
c(Ay) a(x) B(@) Bs + Bi ; 
a. A(Ay) D(a) 
ay = qt {ee Ate PO) ta = ee fh — (10) 


les fonctions propres sont: 
Uy_(4%) = 4 C(An) A(%, An) — %e E(An) 8(% Ax) - 


he Exemples 


a) (x? yl Sey My }(#) ; y(a) oe y(b) = B . (b ily 0) 


Les fonctions propres sont: 


’ V4 Aga 1 x y Z = I! h2 [2- l= Se 
Vy — alle sin are iti = 5 avec 4 ie ae oi , ee In bla ? 
ees ee ee er ee 
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La formule (9) devient ici 


V4 Ay 2 cgay — (—1)k 518 By — Af: (v.-+m2) 


y"+ry = fl*); 
a, (0) + o% 970) =a, Bs vl) + Bay’) = B. (a % SO; Bs By = 0) 


Les fonctions propres sont 


V_ = % SINW, ¥ — A Wy COS MW, * , 


les w, étant les racines positives de l’équation caractéristique 


y wT hy ay W 0. 
Bz cosal— Bywsinwl p,sinal+ By wcosal 
En particularisant (10), nous obtenons la formule 
a On 
a= 47 |nt le — e  - 41| (2x = ©) 
ou 
l af wt— af. 
Na = (Uys Vg) = (02 + 08 wt) = + sin 2 w;,/ — a % sin? @, / 
2 4 ow, 
et 
=e) le Hy we. ig Xe Wh Pe: 
Ck ~ B, cosm, 1 — By, Sinw,!  B, sina, 1 + By wy, Cosa, L 


Remarque. 4 = 0 est valeur propre si les conditions aux limites sont y’(0) = «, 
y’(l) = B, et dans ce cas seulement. Alors Ay) = 0, vy = 1 et 


1 fa-— 
a) — —| os B + Ao]. 


6. Application 4 l’équation de la chaleur 


Pour donner un exemple d’application 4 une équation aux dérivées partielles, 


considérons |’équation 
fap) ou . 


Ox® =" OF,’ 
a laquelle se raméne celle de la chaleur par un choix convenable des unités, la 
condition initiale 


lim U(#, t) =00 (0 < ¥ <1) 
t—> +0 


et les conditions aux limites 


lim U(#,#) = A(t), lim  ~—= Bit). (t > 0) 
x >+0 x—>l—-9 9% 
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Appliquons la transformation de LAPLacE; nous obtenons le probleme 


» WS) =a) - 


uy" 


—su=0, 4u(0,s) = a(s) 
ou 
u(x, 8) = L{U(#, 1}, als) = L{AW}, 0s) = L{ BO}: 
yee Ou i? 0*u 
in ee aa 


En consultant le tableau de la page 152, on trouve que 


=» a : [er 1) = a(s) = ( 1) 2.b(3)| sin(k — 5) F » 


1\ 272 
Mea aCe al |sin(e-4) 2 6, (Oe) 


Nous avons obtenu ce résultat sans recourir a la théorie des résidus. 


7. Conclusion 


Tout probléme aux limites complétement non homogéne est réductible, on 
le sait, A un probléme semi-homogene. L’intérét des séries que nous avons obtenues 
réside dans le fait qu’elles donnent la solution des problemes non homogénes sans 
qu’il soit nécessaire de faire cette réduction, ni de chercher l’intégrale générale de 
l’équation du probléme considéré. Ces séries conviennent parfaitement a l’intégra- 
tion des équations dont le second membre est une fonction génévalement continue, 
équations qui se présentent fréquemment dans les applications. 

Le tableau a la page 152 donne l’expression des coefficients de la série dans 


quelques cas particuliers. 


Summary 


The purpose of this article is not to use a certain transformation to solve 
boundary problems, but to find an operator (viz. a type of series) adapted to the 


given problem. 
In this paper, we limit ourselves to problems for which the homogeneous 


corresponding problem is self-adjoint; the method is applied to problems of the 
second order and to a partial differential equation. 

The obtained series give the solution of completely non-homogeneous pio- 
blems without our having to reduce them to semi-homogeneous problems or to 


find the general solution of their equation. 


(Regu le 27 novembre 1952.) 
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152 


a AG 


pis (4 
a [J 4(1—) — »] eed) ted 


il 1 


(0 + ¥) 


[te + at )— 2] LI. 


[° + ie ta 


oe ee ea ee) (0 + ™) 
(1) 0=70% |g =(1),4 + (Na 
«+ yus T z\ 2 ” = (0),4 *0 + (0)4 To 
i ) by RA a 
vi(Z {) aS aac )| O=-7@399 | J = (),4 ‘n= (0)4 
2/22 ae Scale, =] 
ae es z)s09 iPac )| 0 =1@ 330 g = (4 in = (0),4 
xt ys09 0=10% |g -(),4 m= (0),4 
xt yus 0 = 1 34 g = (1)4 :» = (0)a 
b 
4a ae suol}Ipuos 


(*){ = 44+ ,4 uoyonbz 
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Uber die Instabilitat beim Verfahren der zentralen Differenzen 
fiir Differentialgleichungen zweiter Ordnung 
Von LotHar CoLiatz, Hamburg?) 
Bei SOME Ew ciee numerischer Integration gewohnlicher Differentialgleichun- 
gen y() = f(x, y, y’, ..., y™-1)) kann es passieren, dass ein an sich gutes Verfahren 


mit kleinem Fehlerglied in der numerischen Brauchbarkeit durch ein instabiles 
Verhalten?) («Aufschaukeln», «Aufrauhung») mehr oder weniger stark beeintrach- 


B B 


Al max =40_ 
14 


Fig. 1 


tigt wird. In dieser Zeitschrift hat H. RUuTISHAUSER’) eine Erklarung dieser Er- 
scheinung gegeben; er vergleicht das Anwachsen kleiner Stérungen nach & Schrit- 
ten der Schrittweite #, welches unter der Annahme konstanter partieller Ablei- 
tungen von f nach y, y’,..., yi—}) durch | 4|* gegeben wird, mit dem Anwachsen 
der Lésungen der Differentialgleichung, welches etwa durch |A|* gegeben sei. 
Dabei ist 7 die dem Betrage nach grésste Wurzel einer gewissen algebraischen 
Gleichung. Durch Verwendung genauerer Naherungsformeln kann sich der Grad 
der algebraischen Gleichung erhdhen, und es werden weitere 4 «eingeschleppt», 
fiir welche | 4|>|A| sein kann und welche dadurch Instabilitat bewirken konnen. 

Diskutiert man die Verhaltnisse noch etwas weitergehend, zum Beispiel fur 
das Verfahren der zentralen Differenzen bei Differentialgleichungen zweiter 
Ordnung y” = f(*, y, y’), so ergibt sich das folgende Bild, welches vielleicht von 


1) Forschungsstelle fiir praktische Mathematik der Universitat Hamburg. 

2) Auf diese Erscheinung wurde hingewiesen in L. CoLLarz und R. ZurmUn., Beitrage su dem 
Interpolationsverfahren der numerischen Integration von Differentialgleichungen erster und zweiier 
Ordnung, Z. angew. Math. Mech. 22, 42-55 (1942), speziell Seite 46. 

3) H, RUTISHAUSER, Uber die Instabilitat von Methoden zur Integration gewohnlicher Differential- 


gleichungen, ZAMP 3, 65-74 (1952). 
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Interesse ist, da das Verfahren der zentralen Differenzen wegen mancher Vorziige 
ein vielbenutztes Verfahren ist. 
Bei diesem rechnet man nach den Formeln B 
i 


he 
Veer = 2 Vn — Vea + 2 (fern + 10% pa) 


h 
View = Vea + oor fea ot i ee 


Mit dem Ansatz fiir Stérungen 9, = pp A", np = | Ay hs 
kommt man mit den Abkiirzungen A = h? f, und 
B=hf,, auf die Gleichung (hergeleitet bei Rutis- 
HAUSER; dort ist lediglich a = A/12 und b= B/3 
geschrieben; wir beschranken uns hier auf das Mit- 
teilen der Ergebnisse und verweisen beziiglich der 
Rechnung auf RUTISHAUSER) 


(A~1)[12 (481), 4 1) — ASE 1) OP 1 
—4B(A—1) (424 4441) =D: 


‘Die Werte A berechnen sich mit A= e@ aus 
o?— Boa— A= 0. Die Figuren 1 und 2 geben in 
einer (A, B)-Ebene die Werte des maximalen | /| 
bzw. des maximalen |A|. Figur 3 zeigt in einer 
(A, B)-Ebene zwei Gebiete I, II; in Istimmen maximales |/| und | A| mit grosser 
Genauigkeit iiberein; fiir das Verfahren der zentralen Differenzen ist dort keine 
Instabilitat zu befiirchten; im Gebiet II dagegen besteht Instabilitatsgefahr, so- 
fern man nicht vorbeugende Massnahmen ergreift?). 

Ahnliche Uberlegungen lassen sich auch bei verschiedenen anderen Verfahren 
durchfiihren. 


=i! 
Fig. 3 


Summary 


In the numerical solution of an ordinary differential equation by step by step 
methods it may be that the computational solution is instable and has after a 
certain number of integration steps of the length # nothing to do with the solution 
of the differential equation. For the method of central differences for the diffe- 
rential equation v” = f(x, y, y’) figure 3 gives in the (A, B)-plane (with A = h? Te 
and B= hf,,) a general view of the risk of instability; in the region I there is 
no risk for instability on account of the method of central differences. 


(Eingegangen: 25. September 1952.) 


1) Vgl. zum Beispiel L. Cortatz, Numerische Behandlung von Differentialgleichungen (Springer, 
Berlin 1951), S. 51. 
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Beitrage zur numerischen Behandlung der Schrédinger -Gleichung 
im Falle des Yukawa-Potentials 


Von KALERVO V.LAURIKAINEN und ERKKIK. EvuraAnto, Turku, Finnland?) 


Die Eigenwerte der radialen Schrédinger-Gleichung 


ao [4 1+) ,e 
ae [Fe ee b= 0 (1) 


mit den Randbedingungen ! 
®(0) = B(co) = 0 (2) 


sind neulich im Falle 7 = 0 von HuLTHEN und LAURIKAINEN durch Benutzung 
der Rayleigh-Ritzschen Variationsmethode tabuliert worden [1]}?). Fiir die Figen- 
funktionen wurde dabei der von HULTHEN vorgeschlagene Ansatz 


n 
D =e (1%) Sher (3) 
v=0 


angewandt [2]. Fiir grosse Werte von k wurden die drei niedrigsten Eigenwerte 
b,, bg, bs, unabhangig von der Variationsmethode, auch zu semikonvergenten 
asymptotischen Potenzreihen nach Potenzen von 1 = 1/(2k +2) entwickelt. 
(k wurde als gegeben, 0 als Eigenwertparameter betrachtet.) 

Wir haben versucht, diese Berechnungen fiir andere Werte von / zu verall- 
gemeinern. Zuerst war es moglich [3], verallgemeinerte asymptotische Entwick- 
lungen fiir grosse Werte von & zu gewinnen, wenn l eine beliebige ganze Zahl 
(= 0) ist. Wahlt man 

1 


2k + 2( +1) 


als Entwicklungsparameter, so lautet die Entwicklung des ersten Eigenwertes: 


(4) 


t= 


Ped _ G4 nrc Sens (2iet 4i+ 3) 0 
Tt 


Pian, 
— 5 (+ 1)8 (1018 + 247+ 15) 7° 


(4 1)2 (464 74 + 946 9 + 664 1? + 7611 + 495) a 


1 
15 

4 + (} 1)2 (4560 1 + 23360 1 + 24280 4 — 37020 — 881201 
90 


— 606057 — 15015) 4+::-. 


Wenn man den Variationsansatz (3) fiir positive /-Werte verallgemeinern will, 
muss man vor allem das Verhalten der Eigenlésungen bei grossen und kleinen 
Werten von x beachten. Wenn ¥ -> oo geht, so reduziert sich (1) asymptotisch 


in die Gleichung 
0 _ [ys EM) onc. o 
ax a 


1) Universitat Turku. 


: a IE ic GO SS 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 158. 
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Die im Unendlichen verschwindende Lésung ist hier . 


1 d U+1 . 
= const #/+1 (= : az) gare. (¥ > oo) (7) 


2 dx 


as 
und dasselbe asymptotische Verhalten muss also die Ansatzfunktion aufweisen. 
In der Nahe von ¥ = 0_haben wir entsprechenderweise folgenden asymptotischen 


Ausdruck 


®., = consta™ . (~ > 0) (8) 


as 


Wir haben zunachst den Fall 7 =1 untersucht, wobei insbesondere die fol- 
genden Ansatze versuchsweise gebraucht wurden: 


n 
@ —= eke (1 - ee: a h, e-ve, (9). 
n 
= eke (k + =) (1 — e-®)8 Sh, e* (10) 
vy=0 
(vergleiche zum Beispiel [4], IV C, (119), S. 412), 
L—e? * 
@ = eke (: a —) (Li 62") * oO ieee on (ty 
* v=0 
Man wird erwarten, dass (10) und (11) eine bessere Konvergenz fiir die Variations- 
methode erzeugen als (9), weil erstere sich fiir grosse Werte von x den Eigenfunk- 
tionen naher anschliessen. Fiir kleine Werte von & gibt (10) in der Tat befriedi- 


gendere Resultate als (9), fiir grosse Werte von & verhalt es sich aber gerade um- 
gekehrt, wiaihrend (11) fiir alle k-Werte die besten Resultate liefert (Tabelle 1). 


Tabelle 1 


Vergleich dey Ansatzfunktionen (9), (10) und (17) 
b, berechnet mit je 2, 3 und 4 Parametern h, 


Anzahl 
Ansatz | der Para- 
meter h, 


2 9,186.45, 10,153 06, 11,923 70, 19,79910, 
3 9,186 36, 10,15068, 11,92069, 19,797 02, 
4 9,16103, 10,14973, 11,92051, 19,797 00, 
(10) 2 9,082.83, 10,155 744 11,94789, 19,96377, 
3 9,08209, 10,14927, 11,921 24, 19,814.65, 
4 9,081 98, 10,14912, 11,92046, 19,798 76, 
(11) 2 9,082.83, 10,14977, 11,921 33, 19,798 12, 
3 9,08209, 10,14915, 11,92048, 19,797 01, 
4 9,081 98, 10,14912, 11,92046, 19,797 00, 


Im Falle 7 = 2 wurde dann der Ansatz 


i 1 ma ex 1 ae e7# 2 n 
=k 2 ee ES 
Dm orto hd 4 3h 5 +3 ~ ) | a= es She (12) 
gebraucht. Die Konvergenz des Variationsverfahrens ist dabei sehr befriedigend, 
aber die Berechnungen werden schon ziemlich miihsam. Die Ansatzfunktionen 


/ 
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Ord- 


nungs- | Anzahl 
zahl des | der Para- 
Eigen- | meter h, 


wertes 


bs 


Aw 
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Tabelle 2 
Beispiele dev approximativen Eigenwerte 
p= k = 0,5 Be pas 
a) 1 = 0; Ansatz (3); drei erste Eigenwerte 
1,67993, PKS) 3,81663, 7,896 63, 
1,67985,. ZOO 23> 3,81662, 7,89661, 
1,67981,. 2,769 22, 3,81662, 7,89661, 
6,448 49) Sse, 10,96775, 19,345 21, 
6,447 71; 8,77057; 10,964 04, 193342015 
14,37. Wy Pe, 21,3828, 34,210, 
b) 7 = 1; Ansatz (11); zwei erste Eigenwerte 
9,082 83, 10,149 77, IASOZ1 332 19,79812, 
9,082 09, 10,14915, 11,92048, 19,79701, 
9,081 98, 10,14912, 11,920 46; 19,797 00, 
17,76476, 19,668 52, 22,561 34, 34,797 90. 
17,745 52, 19,65087, 22,93965, 34,775 36 
c) 1 = 2; Ansatz (12); erster Eigenwert 
21,904 61, 22,864 64, 24,968 98, 36,017 63, 
21,895 24, 22 SII 355 24,964.00, 36,011 50, 
21,895 00; 22,859 26, 24,963 93, 36,01147, 


ee ee ee ee ena 


(3), (11), (12) 


immer gentigende Konvergenz errel 


grésseren /-Werten noch beschwerlicher. 
Fiir den Fall 7 = 0 geben wir noch Interpolationsformeln fiir die Berechnung 


des ersten Eigenwertes b, und der entsprechenden Koeffizienten ,, wobei wit im 


Ansatz (3) n= 4 gewahlt haben: 


lassen sich fiir alle Werte von / verallgemeinern, wobei ohne Zweifel 
cht wird, die Berechnungen werden aber bei 


b, = 2,44899 4, + 0,10742 744 4 — 0,00025 385 #2 + 0,00000 9812 # 
— 0,00000 0404 #4, 


Ki = SOU - 0,064.96 60 ¢t — 0,000 26 54 #2 + 0,00000 91 # 


+ 0,00000 04 #4, 


ile 3 1 59 .t= 0,04560 25 ¢ — 0,00025 54 #2 4+ 0,00001 25 # 


+ 0,00000 04 #4, 


eo ODS at 0,005 46 84 ¢ — 0,00000 10 ¢2 — 0,00000 84 2 


+ 0,00000 08 é, 


Hi = MO pee SM SP 0,00590 00 ¢ + 0,00020 46 #2 — 0,00000 09 # 


— 0,00000 46 #4, 


hg = 0;187 1, + 0,004.44 66 ¢ — 0,00014 17 t2 + 0,00000 34 # 


+ 0,00000 17 #, 


(13) 


Ae a aoe 


nee age 
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wo 
k — 0,35 
ur ee are o a4 14 
—— (14) 


ist. Die approximative Eigenfunktion ist hier auf 1 normiert worden. Den Eigen- 
wert b, erhalt man aus (13) im Intervall (0,25... 0,45) mit etwa sieben Dezimalen, 
und in den Parameterwerten h, ist die fiinfte Dezimale mit ein paar Einheiten 
unsicher. 

Eine ausfiihrlichere Darstellung wird etwas spater in den «Annales Universi- 
tatis Turkuensis» ver6ffentlicht [5]. 
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Summary 


Earlier calculations by L. HULTHEN and K.V.LAURIKAINEN concerning 
approximate solutions of the eigenvalue problem (1) and (2) in the case / = 0 [1], 
have been generalized for other values of /. Asymptotic expansions, valid for all 
integer values of /, have been deduced for large values of k [3], the expansion 
for the lowest eigenvalue being (5). As possible generalizations of the function (3)— 
which has turned out to be the ‘‘best fit’? in the case / = 0—for the case / = 1, 
trial functions (9), (10), and (11) have been studied, (11) giving the best results 
(Table 1). In the case /= 2, the trial function (12) was then chosen and was 
found to lead to a very satisfactory convergence. These trial functions can be 
generalized for all values of /, no doubt leading to a good convergence in the 
variational procedure; the calculations, however, become very laborious for 
higher values of /. Some results are given in Table 2. For calculation of the lowest 
eigenvalue b, and the corresponding parameters hf, there are also given inter- 
polation formulas (13) for the case /= 0, taking m = 4 in the trial function (3) 
and normalizing the approximative eigenfunction to unity.-A more detailed 
paper will be published in the Annales Universitatis Turkuensis [5]. 


(Eingegangen: 8. November 1952.) 
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Zur Frage des Dunkelstromes in Photomultipliern 


Von N. ScHAETTI und W. BAUMGARTNER, Ziirich?) 


In einer Arbeit?) tiber den Dunkelstrom von Li-Sb-Photozellen mit Sekundar- 
elektronenvervielfachern konnte gezeigt werden, dass dieser Dunkelstrom keinen 
von der Belichtung unabhangigen Wert besitzt. Eine Komponente des Dunkel- 
stroms, die thermische Emission der Photokathode, ist nach Belichtung des 
Photomultipliers grésser und fallt erst nach einer gewissen Zeit auf den Anfangs- 
wert ab. Inzwischen sind analoge Verhaltnisse bei Vervielfachern mit Cs-Sb- 
Kathoden konstatiert worden. Der Effekt ist hier noch ausgepragter. 


I 
107A 
1,0 
Spannung: 1,5 Volt 
0,9 Messung bei Zimmertemperotur 
x 
\ o--—0 120" 
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08 1\* 
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Fig. 1 
Verlauf der Leitfahigkeit einer Cs-Sb-Schicht in Funktion der Vorbelichtungsdauer. Beginn der 
Messung: 1 min nach Belichtung der Schicht. 


wurden so durchgefiihrt, dass der Vervielfacher bei spannungs- 


Alle Messungen 
m mit einer W-Lampe vorbelichtet und der Dunkelstrom an- 


losem Dynodensyste 
schliessend verfolgt wurde. ped : 
Die Ergebnisse lassen sich wie folgt zusammentiassen . 
il. Durch Belichtung wird die Nullstosszahl einer Li-Sb- und Cs-Sb-Photo- 


kathode erhéht. Die Erhéhung der Nullstosszahl in Funktion der Belichtungs- 
dauer zeigt Sattigung. 


1) Institut fiir technische Physik an der ETH. | 
2 iG soars und W. BAuMGARTNER, Untersuchungen tiber den Dunkelstrom von Photoszellen 


mit Sekunddrelektronenvervielfachern, Helv. Phys. Acta 25, 605 (1952). 
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2. Nach Vorbelichtung der Photokathode findet man eine hdhere Rotempfind- 
lichkeit. Diese klingt parallel der Nullstosszahl im Dunkeln wieder ab. 

3. Das Abklingen der Nullstésse einer vorbelichteten Photokathode kann 
durch Belichtung mit Rotlicht beschleunigt werden (Ausleuchteffekt). 

4. Die spektrale Abhangigkeit von Nullstosszahlerhéhung und Photoeffekt ist 
fiir 400 mp <A < 800myp im wesentlichen gleich. In der Nahe der langwelligen 
Grenze des Photoeffektes aber fallt die Nullstosszahlerhéhung langsamer ab. 

Diese Untersuchungen sind weiter durch Messungen der Leitfahigkeit von 
Cs-Sb-Kathoden nach Vorbelichtung erganzt worden. Diese Kathoden wurden in 
einer zylinderférmigen Zelle zwischen zwei Silberringen formiert1). Wiederum 
wurden sie ohne Spannung vorbelichtet, und anschliessend wurde der Verlauf 
der Schichtleitfahigkeit verfolgt. 

Figur 1 zeigt das Resultat dieser Messungen: Eine Vorbelichtung der Photo- 
kathode erhéht deren Leitfahigkeit. Durch geniigend lange Vorbelichtungsdauer 
wird ein Sattigungswert der Leitfahigkeit erreicht. Die Abhangigkeit dieser Er- 
scheinung von der Lichtwellenlange der Vorbelichtung zeigt einen der Nullstoss- 
zahlerhéhung analogen Verlauf. Insbesondere stimmen die Werte der langwelligen 
Grenze fiir beide Effekte tiberein (~ 775 my fiir Cs—Sb). 

Das obenerwaéhnte Verhalten von Li-Sb- und Cs-Sb-Photokathoden zeigt, 
dass sie neben dem eigentlichen Photoeffekt ebenfalls Eigenschaften aufweisen, 
wie sie bei Photohalbleitern und Kristallphosphoren beobachtet werden. Lage 
und Aussehen des Absorptionsspektrums erinnern speziell an das vielfach unter- 
suchte CdS?). Weitere Untersuchungen sollen zeigen, inwieweit die dort vorge- 
schlagenen Modellvorstellungen beziiglich optischer und elektrischer Eigen- 
schaften auch auf diesen Fall zutreffen. 


Summary 


The investigations on dark current in photomultipliers with Li-Sb- and 
Cs-Sb-photocathodes and on conductivity of Cs-Sb-cathodes indicate, that these 
cathodes show, beside the normal photoelectric effect, the behaviour of some 
photoconductive cells and cristal-phosphors. 


(Eingegangen: 14. Februar 1953.) 
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X. Generalversammlung der Union Radio-Scientifique Internationale 
(U.R.S.I1.) 


Im gegenwartigen Turnus vereinigt sich die U.R.S.I. alle zwei Jahre zu einer 
Generalversammlung. Als erste der elf wissenschaftlichen Unionen der UNESCO 
hat sie sich diesmal nach Australien begeben, wo sie als Gast der australischen 
Commonwealth Scientific and Industrial Research Organisation ihre X. General- 
versammlung durchfiihrte. Der amtierende Prasident, Sir Epwarp APPLETON, 


1) N. Scuarrtr und W. BaumGartner, Cathodes photoélectriques Lithium-Antimoine, Le Vide 
6, 1041 (1951). 


2) J. FAssBENDER, Ann. Phys. [6] 5, 33 (1949). — I. BRoser und R.Warminsky, Ann. Phys. [6], 
7, 289 (1950). 
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leitete die sehr gut besuchte Versammlung — allein aus iiberseeischen Landern 
waren 55 Delegierte eingetroffen, sozusagen alle auf dem Luftwege. 

Die eigentlichen Verhandlungen, wofiir die Universitat Sydney ihre Raum- 
lichkeiten zur Verfiigung stellte, erstreckten sich auf_den Zeitraum vom 12. bis 
21. August. Der Stoff war auf sechs Kommissionen verteilt: 

I. Messtechnik und Normungsfragen, 

II. Wellenausbreitung in der Troposphare, 
III. Wellenausbreitung in der Ionosphare, 
IV. Atmospharische Stérungen, 

V. Radio-Astronomie, 
VI. Schwingungslehre und Schalttechnik. 

Ganz allgemein fiihrten die Verhandlungen zu einer Reihe von Resolutionen 
wissenschaftlicher und organisatorischer Natur. Diese werden demnachst in den 

_Verhandlungsberichten, durch das Generalsekretariat in Brtissel, im einzelnen 

veroffentlicht. Auf den Abdruck der 300 wissenschaftlichen Beitrage — jeweils 
Vol. II der Verhandlungsberichte — wird dagegen verzichtet. Statt dessen sollen 
inskiinftig in vermehrtem Masse Sonderberichte herausgegeben werden. Zu ge- 
gebener Zeit und von Fall zu Fall sollen damit bestimmte, wichtige Fragen in 
ihren weiteren Zusammenhangen tibersichtlich dargestellt werden. 

Den Abschluss der Versammlung bildete ein zweitagiger Besuch in der Bun- 
deshauptstadt Canberra, verbunden mit einem Abstecher nach dem benach- 
barten Commonwealth Observatory auf dem Mt. Stromlo. 

Die folgende, XI. Generalversammlung wird im Jahre 1954 in den Nieder- 
landen durchgefiihrt. Der Prasident fiir die neue Amtsperiode ist R. P. PIERRE 
Leyay (Frankreich). W. GERBER 


Schweizerische Physikalische Gesellschaft 


Die Friihjahrssitzung der Schweizerischen Physikalischen Gesellschaft findet 
am Samstag, dem 2. Mai 1953, im Institut fiir Physik der Universitat Genf statt. 
Anmeldungen fiir Referate aus dem Gebiet der angewandten Mathematik und 
Physik sind umgehend zu richten an die Schweizerische Physikalische Gesell- 
schaft, Sekretariat fiir angewandte Physik (P. pe HALLER, c/o Gebriider Sulzer 
AG., Winterthur), wo weitere Auskiinfte iiber die Veranstaltungen und Besich- 
tigungen wahrend der Tagung zu erhalten sind. Epp HALUER 


4. Kongress der Internationalen Vereinigung fiir Briickenbau und 
Hochbau, Cambridge und London 1952 


Der wissenschaftliche Teil des 4. Kongresses der IVBH. wurde in den Raumen 
der Ingenieurfakultat der Universitat Cambridge durchgefiihrt. Tn sechs Arbeits- 
sitzungen befasste sich der Kongress mit den Themen, deren wichtigste Schluss- 
folgerungen im folgenden zusammengefasst wiedergegeben werden. 


A. Allgemeine Fragen 


I. Bemessungsgrundlagen und Sicherheit. Eine zuverlassige Bemessung eines 
Bauwerkes stiitzt sich sowohl auf die Erfassung aller Belastungen 1m weitesten 
Sinne und auf die méglichst gute Kenntnis der Festigkeitseigenschaften der Bau- 


ZAMP. IV/11 
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stoffe; die dabei noch vorhandenen Liicken sind vor allem durch die Versuchs- 
forschung zu schliessen, wobei eine internationale Zusammenarbeit wiinschbar ist. 

Bei der Behandlung der Frage der Sicherheit spielen Wahrscheinlichkeits- 
betrachtungen eine immer grdssere Rolle. 

Il. Entwicklung der Berechnungsmethoden. Die Berechnungsmethoden passen 
sich dank der Vielfaltigkeit der neueren Untersuchungsverfahren immer besser 
den verschiedenartigsten Problemen an. Die Fortschritte sind besonders augen- 
fallig bei den numerischen Methoden und der Entwicklung besserer experimen- 
teller Verfahren. 


B. Stahlbau 


I. Grundlagen. Im Vordergrund standen hier die Probleme, welche sich aus 
der Anwendung des Schweissens ergeben, wie etwa die Einwirkungen der Tem- 
peratur und der gefahrlichen inneren Spannungen. 

Il. Praktische Anwendungen. Durch die Anwendung von Leichtprofilen, be- 
sonders im Hochbau, durch die geschickte Verwendung von Leichtmetallen oder 
durch geeignete Verfahren zur Vorspannung von Fahrbahndecken bei Verbund- 
tragern sind in der letzten Zeit Fortschritte erzielt worden. Aber auch die klassi- 
schen Konstruktionsarten sind durch zweckmassige Massnahmen zur Aufnahme 
der Krafte noch verbesserungsfahig. Durch eine bessere theoretische und experi- 
mentelle Erfassung der Beulprobleme kénnen die Sicherheitsgrade beziiglich 
Ausbeulen noch weiter verkleinert werden. 


C. Massivbau 


I. Grundlagen und Eigenschaften des Betons. Neben den eigentlichen Festig- 
keitseigenschaften wurden auch allgemeine Betoneigenschaften, wie Regelmassig- 
keit des Gefiiges, Frostbestandigkeit und Verarbeitbarkeit, behandelt. Der Begriff 
der Dichte und seine verschiedenen Bedeutungen erforderte eine Klarstellung. 

Bei den Verwitterungsvorgangen fallt vor allem die Vielfaltigkeit der még- 
lichen Ursachen auf. Von grossem Einfluss kann das Zusammenwirken von sehr 
reinem oder sdurehaltigem Wasser mit Meerwasser oder Meeratmosphare sein. 

II, Aktuelle Probleme des Betons; vorgespannter Beton. Bei der Berechnung 
des Eisenbetons ist die Tendenz festzustellen, die Bruchspannungen an Stelle der 
zulassigen Spannungen bei Gebrauchslast als massgebende Grésse einzufiihren. 

In letzter Zeit wird der vorgespannte Beton auch bei statisch unbestimmten 
Systemen angewendet. Das fiihrt zu Problemen, deren Schwierigkeit besondere 
Anforderungen stellt, die jedoch systematisch einer Lésung entgegengefiihrt 
werden k6nnen. Beim vorgespannten Beton ist bei den Berechnungsmethoden 
und den Konstruktionsverfahren ein Héchstgrad von Genauigkeit und Sorgfalt 
erforderlich, um nicht nur die Tragfahigkeit der Bauwerke, sondern auch deren 
Dauerhaftigkeit zu gewahrleisten. H. WANZENRIED 


Wissenschaftliche Jahrestagung der Gesellschaft fiir angewandte 
Mathematik und Mechanik vom 21. bis 25. April 1953 in Aachen 


Der Vorstandsrat der GaMM hat beschlossen, die diesjahrige Tagung vom 


21. bis 25. April 1953 in Aachen abzuhalten. Anmeldungen sind zu richten an: 
Herrn Prof. Dr. fF. SchuLtz-Grunow, Technische Hochschule, Aachen. 


H. GORTLER 
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Redaktionelles 


Die Redaktionskommission beklagt den Verlust ihres Mitgliedes Herrn Prof. 
Dr. P. Nigcui, der am 13. Januar unerwartet gestorben ist. Die Wiirdigung des 
Lebenswerkes des Verstorbenen wird zu einem spateren Zeitpunkt erfolgen. 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


The Vector Algebra of Vectors and Matrices. By T. L. Wapbr (Addison- 
Wesley Press, Cambridge, Mass., 1951), 189 pp., 10 figs.; $4.50. 


An der Florida State University gebrauchliches Lehrmaterial ist hier zu 
einem Buch von 189 Seiten zusammengefasst, das den Stoff der ersten einer Reihe 
von Vorlesungen tiber Algebra enthalt. 

149 Lehrsdtze und 196 im Text verteilte Ubungsaufgaben kreisen um die 
Grundaxiome der elementaren Vektoralgebra, so dass das Buch besonders fiir 
didaktische Betrachtungen Interesse finden diirfte. Th. Stutz 


Supersonic Flow and Shock Waves. By R. Courant and K. O. FRIED- 
RICHS (Interscience Publishers, New York, 1948) XVI + 464 S., 216 Abb.; $7.00. 


Das vorliegende bedeutungsvolle Buch von R. Courant und K. O. FRIED- 
RICHS will einen Uberblick geben iiber die Uberschallstrémungen und die damit 
verbundenen Fragen, also iiber eines der wichtigsten Gebiete der Aerodynamik. 
Trotz der grossen praktischen Bedeutung dieses Gebietes sind wir heute noch 
weit entfernt von einer allgemein befriedigenden Theorie; ebenso fehlen einheit- 
liche Methoden zur Lésung dér auftretenden Probleme. In den letzten Jahren 
wurden zwar die Bemiihungen sehr intensiviert, vor allem auch in Amerika, 
und manche Fortschritte erreicht, aber das Ziel liegt noch in weiter Ferne. Von 
Bedeutung sind in erster Linie die Untersuchungen der nichtlinearen (quasiline- 
aren) partiellen Differentialgleichungen, die zum Begriff der einfachen Welle 
(simple waves) fiihren, mit deren Hilfe manche schwierigere Str6mungsprobleme 
erfasst werden kénnen. Wie B. RIEMANN schon fand, werden die quasilinearen 
Probleme in der Hodographenebene linear. 

Die Darstellung der Gasdynamik ist in diesem Werke teilweise referierend 
behandelt, indem Probleme nur qualitativ, unter Hervorhebung der Grundge- 
danken, wahrend andere Belange, insbesondere mathematische Entwicklungen 
und Verallgemeinerungen, etwas eingehender ausgefiihrt sind. Fiir Einzelunter- 
suchungen wird auf die einschlagige Literatur (im ganzen sind rund 200 Titel 
aufgefiihrt) hingewiesen, wobei allerdings manche Originalarbeiten in schwer 
zuganglichen Zeitschriften oder in nur mit Miihe erreichbaren Reports veréffent- 
licht sind. Um so dankenswerter ist es, dass R. Courant und K. O, FRIEDRICHS 
versuchten, eine Synthese zu geben, soweit das beim gegenwdrtigen Stand der 
Erkenntnisse iiberhaupt moglich ist. 

Einleitend geben die Autoren eine Zusammenfassung der Thermodynamik 
und der Bewegungsdifferentialgleichungen der kompressiblen Gase und der mathe- 
matischen Theorie der hyperbolischen Differentialgleichungen erster Ordnung 
fiir Funktionen mit zwei und mehr Variablen, und es werden die wichtigsten Be- 
griffe, wie Charakteristik, Hodographentransformation und «simple waves», ein- 
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gefiihrt. Das dritte und langste Kapitel befasst sich mit instationaren, eindimen- 
sionalen Stromungen (Funktionen von *, f) und behandelt Verdiinnungs- und 
Verdichtungsstrémungen, Verdichtungsstésse, sowie in allgemeiner Ubersicht 
das Zusammenwirken mehrerer solcher Erscheinungen. An diesen Entwicklungen 
waren die beiden Verfasser wahrend des Krieges massgebend beteiligt. Anschlies- 
send folgen Untersuchungen von Explosionswellen und Verbrennungsvorgangen. 
Im nachsten Teil kommen die ebenen isentropischen wirbelfreien stationaren 
Strémungen zur Sprache. Als wichtigste hier behandelte Probleme sind zu nen- 
nen: Strémungen durch Kandle und Diisen, Verdichtungsstésse und damit zu- 
sammenhiangende Fragen, Uberschallstromungen um Tragfliigel gemass den St6- 
rungsmethoden. Bemerkenswert ist die hier eingeschaltete allgemeine Diskussion 
iiber das Problem der Randbedingungen bei stationaren Str6mungen. Gas- 
str6mungen durch zwei- und dreidimensionale Kanale und das Ausstr6men 
aus Diisen (speziell Lavaldiisen, wie sie bei Turbinen, Raketen, Windkandlen 
auftreten) werden im fiinften Kapitel behandelt. Das letzte Kapitel befasst sich 
mit dreidimensionalen Str6émungen, insbesondere mit den zylindersymmetrischen 
StrO6mungen um schlanke KGrper, mit stationaren, isentropischen und wirbel- 
freien kegeligen Str6mungen und instationaren Kugelwellen (Explosionswellen 
in Wasser, Luft). 

So gibt das Buch, dessen Text durch zahlreiche Figuren unterstiitzt wird, 
einen vorziiglichen Uberblick iiber dieses interessante physikalische Gebiet, und 
gelegentlich werden auch allgemeine erkenntnistheoretische Fragen gestreift. Die 
Verfasser bemiihten sich, weder zu sehr die mathematischen Probleme noch die 
physikalisch-experimentelle Seite in den Vordergrund zu riicken, so dass gleicher- 
weise Physiker, Ingenieure und Mathematiker dieses Buch mit Gewinn zur Hand 
nehmen werden. E. Roth-Desmeules 


Adhesion and Adhesives. By N. A. DE BRUYNE and R. Houwrnxk (Elsevier 
Publishing Co., Amsterdam, 1951), 518 pp., 205 figs.; 70s. 


Das allgemeine Bestreben der Technik, die Werkstoffe dem jeweiligen Ver- 
wendungszweck méglichst anzupassen, um eine maximale Ausniitzung zu errei- 
chen, bringt es bei der grossen Zahl heute zur Verfiigung stehender Werkstoffe 
mit sich, dass auch Fragen iiber die Verbindungsméglichkeiten der einzelnen 
Werkstoffe miteinander an Bedeutung stark zugenommen haben. Dabei sind vor 
allem bei ungleichartigen Werkstoffpaaren (Metall — Keramik, Metall — Kunst- 
stoffe usw.) die Verbindungen besonders von Interesse, die mit Hilfe eines Binde- 
mittels (Leim, Klebstoff, Kitt, Adhesiv usw.) erreicht werden. Gerade in den letz- 
ten Jahren sind auf diesem Gebiet, wie das Beispiel des bekannten Klebstoffes 
«Araldit» der Ciba zeigt, grosse Erfolge erzielt worden, die erlauben, konstruktiv 
prinzipiell neue Wege zu gehen. Die Technologie des Verbindens stofflich gleich- 
artiger oder ungleichartiger Konstruktionselemente mit Hilfe eines Bindemittels 
ist ausserst komplex, denn es spielen hier vor allem Fragen der Oberflachen- 
beschaffenheit und damit der wirksam werdenden Grenzflachenkrafte eine wesent- 
liche Rolle, die bei der Wahl des Bindemittels und bei der Herstellung der Ver- 
bindung beriicksichtigt werden miissen. Obschon noch viele Einfliisse, welche die 
Haftfestigkeit zweier Grenzflachen bestimmen, nicht in allen Einzelheiten abge- 
klart sind, ist es ausserordentlich zu begriissen, dass in diesem Buch die Gesamt- 
probleme der Haftung auf Grund der in den letzten Jahren erzielten theoretischen 
und praktischen Ergebnissen behandelt werden. Die Verfasser haben bei der 
ausserordentlichen Vielgestaltigkeit des Stoffes den einzig richtigen Weg gewahlt 
und jedes Spezialgebiet durch einen Fachmann behandeln lassen. Was dadurch 
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vielleicht an Kontinuitaét des gesamten Stoffes verlorengeht, wird mehr als wett- 
gemacht durch eine klare und griindliche Darstellung der einzelnen Stoffgebiete. 
Die wichtigsten Kapitel behandeln dabei folgende Themata: 

a) Theoretischer Teil: Die allgemeinen Bedingungen fiir Benetzung von Ober- 
flachen und Adhesionen von Grenzflachen. Molekulare Krafte, die an Grenzflachen 
wirksam sind. Das Verhalten von Klebstoffen wahrend und nach dem Abbinde- 
prozess bei mechanischer Beanspruchung. Kraftlinienverlauf in der Verbindungs- 
schicht bei statischer Beanspruchung in Funktion der Art der Grenzflachen. 
_b) Praktischey Teil: Aufbau und Verwendung organischer Bindemittel, auf- 
geteilt in tierische, pflanzliche, synthetische und asphaltartige Stoffe. Anorga- 
nische Bindemittel, umfassend Zemente und metallische Verbindungen (Lote). 
Spezielle Bindemittel fiir natiirliche und synthetische Kautschuke. Physikalische 
Priifung der Bindemittel und der Haftfestigkeit. 

Am Schluss jedes Kapitels folgt eine ausfiihrliche Literaturiibersicht, wobei 
auch die wichtigsten Patente angefiihrt sind. 

Trotz der Vielgestaltigkeit der in den einzelnen Kapiteln behandelten Probleme 
ist es den Autoren gelungen, eine geschlossene Darstellung des Stoffes zu geben, 
wobei vor allem die Tendenz, allgemein giiltige Zusammenhange hervorzuheben, 
erfreulich stark hervortritt. Die meisten Kapitel sind nicht nur fiir den Chemiker 
von Interesse, sondern ebensosehr fiir den konstruktiv tatigen Ingenieur, sei 
es im Holz- oder Maschinenbau oder in der Elektrotechnik. Uberall stellen sich 
Probleme iiber zweckmassige Verbindung von Werkstoffen, fiir deren Lésung 
das Studium des vorliegenden Buches sicher wertvolle Anregungen geben wird. 
Beziiglich Ausstattung befriedigt der bekannte Verlag auch verwohnteste An- 
spriiche, und das Buch darf als Standardwerk auf diesem speziellen Gebiet der 
Werkstoffe allen interessierten Fachleuten bestens empfohlen werden. fF, Held 


The Magnetron. By R. Latruay, A. H. Krtnc and L. RusHrortH (Chapman 
& Hall, London, 1952). 142 pp., 82 figs.; 18s. 

Die Verfasser sind bekannt durch Beitrage zur Entwicklung des Multi-Cavity- 
Magnetrons in Blockkonstruktion. Das Buch gibt mit einer historischen Einlei- 
tung, welche auch andere Oszillatoren umfasst, eine zusammenfassende Darstel- 
lung des Magnetrons speziell im Hinblick auf die Verwendung fiir Radar auf 
Zentimeterwellen mit Impulstechnik hoher Leistungen. Die allgemeinen Eigen- 
schaften des Multi-Cavity-Systems (gekoppelte Resonanzkreise) und die Lei- 
stungsauskopplung werden ausfiihrlich in zwei weiteren ‘Kapiteln behandelt. 
Drei folgende Kapitel diskutieren den Elektronenmechanismus, Schwingungs- 
unterhaltung (mode stability), Schwingungsanfachung, Wirkungsgrad, ferner die 
Probleme der Kathodenkonstruktion fiir hohe Emissionsstr6me, die Konstruk- 
tion und Fabrikation des Magnetrons. Anschliessend folgt das Kapitel tiber die 
Belastung des Magnetrons mit Betriebscharakteristiken (Belastungsdiagramme, 
Leistungsmessung). Den Schluss bildet die Anwendung des Magnetrons fiir Radar 
mit zugehériger Modulationstechnik. 

Den Verfassern gelingt im ganzen eine anschauliche Darstellung, welche durch 
die 82 schematischen Figuren und Photographien vorteilhaft unterstutzt wird. 
Das Literaturverzeichnis nach jedem Kapitel (hauptsachlich Reports) zeugt von 
grosser Verarbeitung des Materials. Ob den, Verfassern auch die anschauliche 
Darstellung des Elektronenmechanismus (nach den schwierigen Originalarbeiten) 
gelungen ist, bleibe dem Leser zur Beurteilung tiberlassen. Jedenfalls wendet 
sich das Buch nicht nur an den ‘Aussenstehenden, sondern auch der Kachmann 


kann vieles daraus lernen. F. Lid 
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Partielle Differentialgleichungen. Von J. Horn (W. de Gruyter, Berlin 
1949). 228 S., 8 Abb.; DM 14.-. 


Die vorliegende vierte Auflage dieses bewahrten Lehrbuches ist ein Neudruck 
und entspricht damit im wesentlichen der zweiten Auflage aus dem Jahre 1929. 
Einzig das Literaturverzeichnis, das nur das Wichtigste umfasst, ist bis etwa 1940 
nachgefiihrt. Das Buch bietet den klassischen Bestand der Theorie der partiellen 
Differentialgleichungen unter Beschrankung aut die am haufigsten auftretenden 
partiellen Differentialgleichungen erster und zweiter Ordnung mit zwei unabhan- 
gigen Variablen. Die allgemeinen theoretischen Entwicklungen werden klar und 
ausfiihrlich dargestellt, die Existenz der Lésungen wird bewiesen. Meist sind die 
verschiedenen Lésungsmethoden noch an wichtigen Beispielen aus der theoreti- 
schen Physik (Mechanik, Thermodynamik, Potentialtheorie) und der Flachen- 
theorie vollstandig durchgefiihrt. Inhaltsiibersicht: Normalformen der linearen 
partiellen Differentialgleichungen (DGIn.) zweiter Ordnung mit zwei unabhan- 
gigen Variablen; hyperbolische DGln.; Theorie der linearen Integralgleichungen 
nach D. HIrBerT und E- Scumipt, Lésung von FREDHOLM; Randwertaufgaben 
fiir gewohnliche lineare DGln., soweit diese fiir die partiellen DGIn. von Bedeu- 
tung sind; Randwertaufgaben fiir elliptische DGln.; parabolische DGln.; par- 
tielle DGln. erster Ordnung mit zwei unabhangigen Variablen; partielle DGln. 
zweiter Ordnung mit zwei unabhangigen Variablen (Monge-Ampéresche DGl.). 

Das gut und iibersichtlich aufgebaute Werk kann den Studenten der Mathe- 
matik und der Physik sowie der Ingenieurwissenschaften sehr empfohlen werden. 
Allerdings sei abschliessend noch darauf hingewiesen, dass das Buch — und dies 
lag in der Absicht des Verfassers, der nur die «reine» Theorie darstellen wollte — 
nichts enthalt tiber numerische Methoden zur Lésung von nicht geschlossen 
integrierbaren Fallen, die ja in vielen Gebieten eine Rolle spielen und von grosser 
praktischer Bedeutung sind. E. Roth-Desmeules 


Tables of the Error Function and of the First Twenty Derivatives. 
By the Staff of the Computation Laboratory (The Annals of the Computation 
Laboratory of the Harvard University, Vol. 23). (Harvard University Press, 
Cambridge, Mass., 1952). 300 pp.; $8.00. 


Es handelt sich um eine ausfiihrliche, sechsstellige Tafel des Fehlerintegrals 
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und dessen 21 ersten Ableitungen [q()(%) bis g(2(x)], die jedem, der mit dem 
Fehlerintegral arbeiten muss, eine willkommene Hilfe sein wird. 

Die Tabelle reicht so weit, als die Funktionswerte noch die Schwelle von 
10~* iibersteigen, und das Argumentintervall (meist 0,002 oder 0,004) ist so ge- 
wahlt, dass mit Hilfe der mit tabellierten Ableitungen leicht nach HERMITE oder 
auch mit einer Taylor-Reihe interpoliert werden kann. 

Den eigentlichen Tabellen (276 Seiten) geht eine Einleitung voraus, in der die 
wichtigsten funktionentheoretischen Eigenschaften der tabellierten Funktionen 
und der damit zusammenhangenden Hermiteschen Polynome sowie einige An- 
wendungsméglichkeiten mit den zugehérigen Literaturangaben zusammengestellt 
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Geratebauanstalt Balzers 


Zweistufige rotierende 
Vakuumpumpen 


Saugleistungen von 5 bis 180 m® pro Stunde. 
Selbsttatige Olregeneration, einstellbares Gas- 
spulventil zur direkten Absaugung kondensier- 
barer Gase und gesattigter Dampfe, zwangs- 
laufige Luftkiihlung. Kleiner Olinhalt. End- 
vakuum garantiert besser als 2 x 10-4 mm Hg 


Verkaufsorganisation 


Trih, Tauber & Co. AG., Ziirich 


Fabrik elektrischer Messinstrumente und wissenschaftlicher Apparate 


Muba Halle IIIb, 2. Stock, Stand 2640 


Praxis der Physik, Chemie und Photographie 


im Unterricht der Schulen 


Erscheinungsweise monatlich, Umfang 40 Seiten, reich illustriert. 
Bezugspreis im Quartal (3 Hefte) DM 4.50 plus Zustellgebiihr — Fur 
Studienreferendare und Studierende ermassigter Bezugspreis DM 3.15 
je Quartal bei direktem Bezug vom Verlag — Probenummern stehen 
auf Anforderung zur Verfiigung — Im Jahre 1952 wurden unter an- 
derem folgende Original-Veréffentlichungen herausgebracht: 


Physik 
Elektrische Wechselspannungen mittels Staubfiguren — Ultraschall- 
versuche (Katakaustika) — Schulversuche zur Radioaktivitat — Bau 
einer Nebelkammer — lonenwanderung im Schulversuch — Herstellung 
von Polarisationsfiltern im Unterricht — Vektographen (Polarisations- 


anaglyphen) Cara 
i 


Elektronentheorie im Chemieunterricht — Kunststoffe im Unterricht — 
Farbstoffe im Schulversuch — Schaumstoffe — Chemische Unter- 
richtsversuche aus Aufbauteilen — Strukturformel des Athanols 


AULIS-VERLAG Frankenberg an der Eder (West-Deutschland) 
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Werke der reinen und angewandten Mathematik und Physik 


Mathematik 
B1eBERBACH, L.: Theorie der geometrischen Konstruktionen. . . geb. Fr. 18.70 
CarRATHEODORY, C.: Funktionentheorie in zwei Banden Band I geb. Fr. 39.50 


Band II geb. Fr. 27.05 
Caratutopory, C.: Mass und Integral und ihre Algebraisierung, erscheint 1953. 
Doetscu, G.: Handbuch der Laplace-Transformation. Band I: Die theoretischen 


Grundlagen der Laplace-Transformation ..... . ¢ sos geboFrgss 20 
FINSsLER, P.: Kurven und Flachen in allgemeinen Raumen . . .  geb. Fr. 15.40 
Fueter, R.: Analytische Geometrie der Ebene und des Raumes . geb. Fr. 24.95 
Herrmann, H.: Ubungen zur projektiven Geometrie . .... . geb. Fr. 17.— 
Hurwitz,A.: Mathematische Werke, 2 Bande ....... . geb.je Fr. 54.10 
LinpDER,A.: Statistische Methoden, 2. Auflage 1951. ..... =. geb. Fr. 31.20 
LocuEer-Ernst,L.: Differential- und Integralrechnung . . . . geb. Fr. 49.90 


Locuer-Ernst,L.: Einfiihrung in die freie Geometrie ebener Kueven geb. Fr. 12.50 
Osrrowsk1,A.: Vorlesungen tiber Diff.- u. Integralrechng., BandI_ geb. Fr. 36.40 
Band II geb. Fr. 69.70 
ScuLArii,L.: Gesammelte mathematische Abhandlungen, BandI  geb. Fr. 56.15 
BandII = geb. Fr. 56.15 


STIEFEL, E.: Lehrbuch der darstellenden Geometrie. .. . geb. Fr. 29.65 
VoELKER, D., Dortscn, G.: Die zweidimensionale Unpladestvanct: geb. Fr. 47.85 
ZWINGGI, E.: "Versicherungsmathematik ss oe) ole te ccaeie. Oem OD erai te Reene 


Phyotk und Astronomie 


GomsAs, P.: Theorie und peer e eery: des Mehrteilchenproblems der Wellen- 


mechanik .. 2 a SO PebS ieaiciae 
NIETHAMMER, TH.: Die veniaaiak Methoden der asteonomisch: -geographischen Orts- 
bestinamaung: “ts . SMa. chgeb. Pessese 


VON DER PAHLEN, E.: BintViramed in ndie ‘Dynamik ve von Stemware geb. Fr. 39.50 
WALDMEIER, M.: ‘Bintibrune ¢ in die’ Astrophysik. .y\) i.) sa eeb. Reeagae 
WALDMEIER, M.: Tabellen zur heliographischen Ortsbestimmung ._ geb. Fr. 14.55 


WALDMEIER, M.: Die Sonnenkorona, BandI... . ose sw Reb. Fs Sava 

Kongressbericht der Internationalen Konferenz tiber Keraphysik und theoretische 

Physik, Basel 1949, 248 S. mit vielen Illustrationen. Broschiert Fr. 16.65 
Ingenteurwissenscha/len 

Bopea,E.: Giorgis rationales MKS-MaBsystem. Broschiert. . . . Fr, 28.10 

GOLDSTEIN, J.: Die Messwandler. 45 2). <.0 > a ee ee ee Peter soale 

JAEGER, CH: Technische Hydraulik. . ... . . geb. Fr; 52.— 


KuHLMANN, K.: Theoretische Elektrotechnik. Band ll: Geanasiice der Theorie der 
Wechaelstrime und des einphasigen Transformators. . . . . . geb. Fr. 27.55 


Band III: Grundztige der Theorie elektrischer Maschinen. . .  geb. Fr. 74.90 
MetssneEr, E., und Zrecrer, H.: Mechanik. Band I: Statik der starren, fliissigen und 
elastischen Worperimee aves abi: yao ee RN esl CDi Camere 
Band II: Dynamik der Roxen Korper ee ey We ht Oe gh ore hy CRG 
Band III: Dynamik der Systeme... . geb. Fr. 46.80 
MicHaEL, W.: Ortskurvengeometrie in der komplexen Ziiiensbens geb. Fr. 13.50 
SANGER, R.: Ballistische Stérungstheorie . . we hie ot) oo ae tees. 6 CED. BRE LOneD 
Srusst, R.: Vorlesungen tiber Baustatik, Band L of) a WyreeteneG os, oe) BED Pita Ose 
RusTERHOLz,A.: Elektronenoptik, Bandi: Grundziige . . . . . .  geb. Fr. 31.20 


RicuTer, R.: Elektrische Maschinen, Band I (Band II 1953) . . . geb. Fr. 49.40 
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ARCHIV DER MATHEMATIK 


ARCHIVES OF MATHEMATICS ARCHIVES DE MATHEMATIQUES 


Herausgegeben in Verbindung 
mit dem Mathematischen Forschungsinstitut in Oberwolfach 


von H. KNrEsER und W. Stss 


Beirat: G. Bot, E. Bomp1ani, P. TEN BruGGENcATE, J. DiruDONNE, Cu. ExnRES- 
MANN, H. GOrTLER, H. Hapwicer, H. Horr, W. Macnus, T. NaGeLi, Cur. PAvc, 


- J. Ravon, K. ReipeMersTER, J. A. SCHOUTEN, H. SEIFERT, E, SPERNER, E. STIEFEL 
Redaktion: H. BiLHarz, Universitat Wurzburg 


Erscheint seit 1948; ab Band 3 (1952) jahrlich regelmassig in 6 Heften von je ca. 
80 Seiten. Format 17X24 cm. Abonnementspreis pro Band sFr. 60.- (DM 60.-), 
Einzelheft sFr. 12.- (DM 12.-) 


Das ARCHIV DER MATHEMATIK veréffentlicht in erster Linie Originalarbetten aus dem 
Gesamtgebiet der Mathematik und ihrer unmittelbaren Anwendungem (bis zu 10 
j Druckseiten Umfang) in deutscher, englischer, franzdsischer oder italieniseher 
Sprache. Daneben kénnen in beschranktem Masse Selbsireferate ber bislang unver- 

dffentlichte gréssere Arbeiten, deren wissenschaftliche Bedeutung dies gerechtfertigt 

MY : erscheinen lasst, Aufnahme finden. In diesen Selbstreferaten werden ausser den Re- 
: sultaten die wesentlichen Schritte der Beweisfiihrung mitgeteilt. Schliesslich ge- 
langen in zwangloser Folge zusammenfassende Berichte iiber die Fortschritte einzelner 
Sondergebiete, die in rascher Entwicklung begriffen sind, mit ausfihrlichen Litera- 

turangaben zum Abdruck. 
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